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1.  INTRODUCTION 

The  funds  received  from  the  DURP  proposal  were  to  be  used  to  purchase  a  Polytec  scanning 
vibrometer.  The  vibroraeter  is  being  used  to  make  precise  measurements  of  the  mistuned  modes  and 
natural  frequencies  of  integrally  bladed  disks  (EBRs)  in  the  Mistuning  System  Identification  Laboratoiy. 
The  resulting  data  is  then  used  as  input  to  system  identification  software  that  is  being  developed  at  CMU 
under  a  separate  Air  Force  contract.  The  system  identification  software  then  determines  the  mistuning  in 
each  blade.  This  report  documents  the  development  of  the  vibration  laboratory,  the  types  of 
measurements  that  are  being  done,  and  the  research  that  has  been  made  possible  through  its  use. 


2.  EQUIPMENT  PURCHASED  AND  LABORATORY  ESTABLISHED 
2.1  The  Poly  tec  Vibrometer 

The  Polytec  scanning  vibrometer  was  purchased  in  2002.  The  equipment  is  shown  in  Figure  1. 


Figure  1  Pol)^ec  Scanning  Vibrometer 


Some  of  the  important  features  of  the  vibrometer  are: 

•  The  Measurements  Can  Be  Automated 

-  Built-in  Function  Generator  drives  Excitation  Speaker 

-  Automatically  Repeats  Excitation  and  Scans  IBR  for  Modal  Measurements  (several 
hundred  transfer  functions  may  be  required  per  IBR). 

-  More  Precise  Positioning  of  Measurements 

•  Data  Validation  Ensures  High  Quality  Measurements 

-  Automatically  Retakes  Data  if  Signal  is  Not  Optimal 

•  12,800  FFT  Lines 
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-  High  Resolution  Measurements  Needed  Because  of  High  Modal  Density  (e.g.may  have 
50  modes  between  1300  -1350  Hz) 

It  was  important  to  be  able  to  automate  the  measurements  because  we  could  set  up  the  test  and  let 
the  experiment  run  overnight.  Alternatively,  they  would  have  taken  weeks,  if  not  months,  to  run  the 
same  experiment  if  we  had  used  our  previous  equipment. 

Another  advantage  of  the  Polytec  vibrometer  is  its  high  frequency  resolution,  i.e.  12,800  FFT  lines. 
Bladed  disks  often  have  frequencies  that  are  very  close  together  and  a  high  frequency  resolution 
capability  is  required  to  determine  their  frequency  response.  An  example  of  a  representative  frequency 
response  plot  for  an  IBR  is  shown  in  Figure  2.  Note  that  about  ten  modes  are  located  within  a  ten  Hertz 
frequency  band.  Consequently,  the  very  fine  frequency  resolution  capability  of  the  Polytec  vibrometer 
was  needed  to  make  these  measurements. 
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Figure  2  Frequency  Response  of  an  IBR 

A  list  of  the  equipment  purchased  with  DURIP  funds  is  given  in  Attachment  1. 

2.2  Additional  Laboratory  Components 

In  order  to  establish  the  System  Identification  Laboratory  we  have  developed  and  purchased 
additional  equipment  from  other  funding  sources.  The  experiments  are  performed  on  a  vibration  table 
with  base  isolation.  Figure  3  (a).  We  have  a  signal  generation  system  and  amplification  system  that  can 
either  excite  the  blades  acoustically  or  magnetically,  see  Figure  3  (b)  and  Figure  4.  The  excitation 
signals  can  be  phased  from  blade  to  blade  so  as  form  a  traveling  wave  excitation  that  simulates 
excitations  found  in  the  engine,  i.e.  in  the  lab  the  excitation  source  rotates  relative  to  the  IBR  whereas  in 
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the  engine  the  IBR  rotates  relative  to  the  excitation  source.  In  addition,  the  excitations  can  simulate 
transients  and  pulse  excitations.  As  a  result,  we  can: 

•  Run  automated  modal  tests 

•  Identify  the  mistuning  in  the  system  model  the  system  and  predict  the  response  to  a 
traveling  wave  (engine  order)  excitation 

•  Verify  the  accuracy  of  the  predictions  with  a  traveling  wave  experiment. 

This  capability  is  used  extensively  in  the  paper  by  Rossi  et  al.  given  in  the  attachment. 


(a)  Vibration  Isolation  Table  &  Vibrometer  System  (b)  Excitation  Control  System 
Figure  3  Mistuning  System  Identification  Laboratory  Equipment 


Figure  4  Academic  IBR  with  Magnetic  Exciters 
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3.  RESULTING  RESEARCH 

We  have  used  the  Mistuning  System  Identification  Laboratory  extensively  in  our  research  at 
Carnegie  Mellon  University.  It  has  contributed  to  the  research  discussed  in  the  following  papers. 

1.  Feiner,  D.M.  and  Griffin,  J.H.,  2004,  “Mistuning  Identification  of  Bladed  Disks  Using  a 

Fundamental  Mistuning  Model  -  Part  I:  Theory  and  Part  11;  Application”,  Journal  of 

Turbomachinery,  Vol.  126,  pp.  150  -  158  and  159-165. 

2.  D.M.  Feiner,  J.H.  Griffin,  K.W.  Jones,  J.A.  Kenyon,  O.  Mehmed  and  A.P.  Kurkov,  2003, 
“System  Identification  of  Mistuned  Bladed  Disks  from  Traveling  Wave  Response 
Measurements,”  Proceedings  of  the  2003  ASME  Design  Engineering  Technical  Conference, 
Chicago,  IL,  ASME  Paper  DETC  2003AaB-48448. 

3.  Ayers  J.P.,  Feiner  D.M,  and  Griffin,  J.H.,  2004,  “A  Reduced  Order  Model  For  Transient 
Analysis  Of  Bladed  Disk  Forced  Response,”  9th  National  Turbine  Engine  High  Cycle  Fatigue 
(HCF)  Conference,  Pinehurst  ,NC. 

4.  Rossi,  M.R.,  Feiner,  D.M.,  and  Griffin,  J.H.,  2004,  “Experimental  Study  Of  The  Fundamental 
Mistuning  Model  For  Probabilistic  Analysis,”  9th  National  Turbine  Engine  High  Cycle  Fatigue 
(HCF)  Conference,  Pinehurst  ,NC. 

Copies  of  these  papers  are  provided  in  Attachment  2. 
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ATTACHMENT  1 


EQUIIPMENT  PURCHASED  USING  DURIP  FUNDS 


PAGE -6  • 


CAKNEGIE  MELLON  UNIVERSITY 
SHIP  MECHANICAL  ENGINEERING  DEPT. 
TO-  REF:  PO#  127603 

HAMERSCHLAG  HALL  C124 
PITTSBURGH,  PA  15213 

CARNEGIE  MELLON  UNIVERSITY 
TO  ATTN  :  ACCOUNTS  PAYABLE 
5000  FORBES  AVENUE 
PITTSBURGH,  PA  15213-3890 


ORDER  INFORMATION: 
23  Midstate  Drive 
Suite  212 

Auburn,  MA  01501 
PH:  (508)832-3456 
FX:  (508)832-0506 


Invoice  No: 
Shipping  Date: 
Order  Number: 

Oust.  PO#  Date: 
Customer  PO: 

Ship  Via: 

Invoice  Terms: 
Customer  No: 

Sales  Rep: 


REMIT  PAYM^T  TO: 

1342  BELliAVa. 

SUITE  3-a\  1  ^ 

TUSTIN.cA  92f80\A\  0 
PH:  (714)8A5-i]835\ 

FX:  (714)85M8g^\ 


10591 
5/31/02 
127603 
UPS  GROUND 
NET  30  DAYS 
CA0218 
DEO 


QUANTITY 

SHIPPED 

QUANTITY 

.::Elo:: 

/^x  -^MODEL  ' . ■ 
NUMBER  x 

.  DESCRIPTION  . 

'  ■  '  "  '  UNIT  . 

'  ^  ‘  :  price:  :  • 

4-...:- ■ . : . :: 

_ 1 _ 

TQTAt; 

QUANTITY 

ORDERED. 


1  ) 
1  I 
1  1 
1  \ 

1  ^ 
1  I 
1  \ 
1  i 
1  ^ 
1  \ 
1  \ 
1  \ 
1  < 
1  \ 
1  ^ 
1  » 
1 

1  \ 
1  ' 
1  I 
1  \ 
1  \ 


P  PSV-300-H  VIBRASCAN  LASER  VIBROMETE 
0  PSV-300-H-B  PSV-300  HARDWARE  ID-31960 

Q  PSV-300-H-S  PSV-300  SOFTWARE  ID-31981 

0  PSV-Z-082  VIS. BASIC  ENGINE  ID-32548 
O  PSV-Z-061  UFF  SOFTWARE  ID-31057 
O  PSV-Z-062  PROF. GRID  LAYOUT  ID-31616 
P  PSV-SOFTDESK  DESKTOP  VERSION  ID-32609 
O  PSV-Z-066  HIGH  RES.  SCAN  ID-31058 
O  PSV-Z-070-H  ZOOM  FFT  SOFTWRE  ID-32142 
O  PSV-Z-081  HIGH  RES.  FFT  ID-32122J 
O  PSV-i08  HD  TRIPOD  W/MOTR. PAN/TILT 

O  OFV-056-CF99  CLOSE-UP  MODULE  ID-32205 
O  PSV-Z-035  CART  FOR  PSV-300  ID-30639 
O  PSV-Z-018  VERT.  TEST  STAND  ID-32336 
O  PSV-TTT-1  1  DAY  TRAINING/INSTALLATI 
0  PSV-300-SM12  12  MTH.SFT.MAINT. ID-31991 
O  PSV-EW-300-H  1-YR  EXT. HARDWARE  WARRANT 
O  OMB-PSV-4  SET  OF  FOUR  CASES  v/ 

0  PSV-Z-061-D  DSKTOP  VER.W/UFF  ID-3260C> 
O  OFV-CL-80  80MM  CLOSEUP  LNS  ID-31551 
0  OFV-CL-150  SM.  PART  LENS  ID-32226 
Q  OMB-LUZ  RING  LIGHT  ILLUMINATOR 


PSV-300  HARDWARE  CONSIST  OF: 

OFV-3001-SH6  VIBROMETER  CONTROLLER  SN: 5021175 
OFV-056  SCANNING  HEAD  SN: 6020966 
0FV-056-CF99  CLOSE-UP  UNIT  SN: 60120380008 
PSV-PC-H  COMPUTER  CONSOLE  SN: 5021121 
PSV-Z-040-H  JUNCTION  BOX  SN: 50207500004 
PSV-Z-051  REMOTE  FOCUS  HANDSET  SN: 50151800045 
OF^7-303S  SENSOR  HEAD  SN: 6020922 _ ^ _ 


CARNEGIE  MELLON  UNIVERSITY 
SHIP  MECHANICAL  ENGINEERING  DEPT. 
TO:  REF:  PO#  127603 

HAMERSCHLAG  HALL  C124 
PITTSBURGH,  PA  15213 

CARNEGIE  MELLON  UNI^/ERSITY 
TO  ATTN  :  ACCOUNTS  PAYABLE 
5000  FORBES  AVENUE 
PITTSBIJRGH,  PA  15213-3890 


QUANTITY 

QUANTtTY 

QUANTITY 

:  (  MODEL  . 

ordered: 

;r  SHIPPED 

ORDER  INFORMATION: 
23  Midstate  Drive 
Suite  212 

Auburn,  MA  01501 
PH:  (508)832-3456 
FX:  (508)832-0506 

invoice  No: 
Shipping  Date: 
Order  Number: 

Gust.  PO#  Date: 
Customer  PO: 

Ship  Via: 

Invoice  Terms: 
Customer  No: 

Sales  Rep: 


REMIT  PAYMENT  TO: 
1342  BELLAVE. 
•SUITE  3-A 
TUSTIN,  CA  92780 
PH:  (714)850-1835 
FX:  (714)850-1831 

10591 
5/31/02 
127603 
UPS  GROUND 
NET  30  DAYS 
CA0218 
DEO 


bESQRIPTION 


UNIT 
PRICE  . 


;/''  -'.:.:'.-'.;.TOTAL , 


OMB-LUZ  FIBER-OPTIC  RING  LIGH  SN: 126534 
PSV-Z-035  PSV  MOBILE  CABINET  SN;N/A 
VIEWSONIC  VG150  COMPUTER  MINITOR  SN:GQ14851965 
CHERRY  D91275  COMPUTER  KEYBOARD  SN:G026432 
MICROSOFT  INTELLIMOUSE  COMPUTER  MOUSE  SN: 5354164 
PSV-300-H  PSV-300  HARDLOCK  SN:  RUS-00000219 
PSV-DESKTOP  DESKTOP  HARDLOCK  SN:RUS-0000021A 
PT570-24P  MOTORIZED  PAN/TILT  HEAD  SN: 01036-14-0012 
PSV-108  HEAVY  DUTY  TRIPOD  SN:132X/MB28 
PSV-Z-017  PAN  /TILT  HEAD  INTERFACE  PLATE  SN:N/A 
PSV-Z-020-A  TRIPOD  ADAPTER  PLATE  SN:N/A 
PSV-Z-018  VERTICAL  TEST  STAND  SN: 60209480002 
PSV-Z-018(C0NT'D)  CONTROL  UNIT  SN: 3255312 
OMB-PSV-4  SET  OF  4  PSV  CASES  SN:N/A 


Comment: 

GARY  HAYDEN  412-268-6248  ,/  362--17M 


COMMENT  ON  COSTS 


There  were  two  complications  that  initially  affected  the  cost  of  the  vibrometer. 

1.  We  received  several  different  invoiees  from  the  manufacturer  of  the  Polytec  vibrometer. 
Initially,  they  included  charges  for  shipping  that  were  not  consistent  with- their  original 
quote  and  that  caused  the  cost  to  exceed  the  funds  provided  by  the  DURIP  grant. 
Eventually,  we  worked  out  the  problem  so  that  the  actual  cost  matched  their  quote. 

2.  CMU  normally  charges  overhead  on  maintenance  contracts  and  the  purchase  of  the  Polytec 
vibrometer  included  $8,285  for  a  maintenance  contract.  The  overhead  charge  was  $4,143. 
The  cost  of  the  CMU  overhead  was  not  included  in  the  DURIP  proposal  and  this  would 
have  cause  a  cost  over  run.  However,  it  was  pointed  out  to  the  CMU  finance  person  that  the 
cost  of  the  maintenance  contract  was  included  in  the  original  proposal  and  that  CMU  had 
signed  off  o  n  the  b  udget  w  ithout  a  n  o  verhead  c  harge.  A  s  a  r  esult,  C  MU  a  greed  t  o  c  ost 
sharing  the  indirect  (overhead)  costs. 

A  copy  of  the  “Project  Detail  Report”  for  January  2004  is  provided  on  the  next  page  and  documents 
that  the  total  cost  was  equal  to  the  grant  amount  of  $213,630. 
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ATTACHMENT  2 


Relevant  Publications 


D.  M.  Feiner 
J.  H.  Griffin 


Mistuning  Identification  of  Biaded 
Disks  Using  a  Fundamental 
Mistuning  Model— Part  I:  Theory 


Department  of  Mechanical  Engineering, 
Carnegie  Mellon  University, 
5000  Forbes  Avenue, 
Pittsburgh,  PA  1521? 


This  paper  is  the  first  in  a  two-part  study  of  identifying  mistuning  in  biaded  disks.  It 
develops  a  new  method  of  mistuning  identification  based  on  measurements  of  the  vibra¬ 
tory  response  of  the  system  as  a  whole.  As  a  system-based  method,  this  approach  is 
particularly  suited  to  integrally  biaded  rotors,  whose  blades  cannot  be  removed  for  indi¬ 
vidual  measurements.  The  method  is  based  on  a  recently  developed  reduced  order  model 
of  mistuning  called  the  fundamental  mistuning  rfiodel  (FMM)  and  is  applicable  to  isolated 
families  of  modes.  Two  versions  of  FMM  system  identification  are  presented:  a  basic 
Version  that  requires  some  prior  knowledge  of  the  system's  properties,  and  a  somewhat 
more  complex  version  that  determines  the  mistuning  completely  from  experimental 
data.  [DOI:  10.1115/1.1643913] 


1  Introduction 

Biaded  disks  used  in  turbine  engines  are  nominally  designed  to 
be  cyclically  symmetric.  If  this  were  the  case,  then  all  blades 
would  respond  with  the  same  amplitude  when  excited  by  a  trav* 
eling  wave.  However,  in  practice,  the  resonant  amplitudes  of  the 
blades  are  very  sensitive  to  small  changes  in  their  properties. 
Therefore,  the  small  variations  that  result  from  the  manufacturing 
process  and  wear  cause  some  blades  to  have  a  significantly  higher 
response  and  may  cause  them  to  fail  from  high  cycle  fatigue.  This 
phenomenon  is  referred  to  as  the  mistuning  problem,  and  has  been 
studied  extensively.  Srinivasan  provides  a  thorough  review  of  this 
topic  in  [1]. 

To  address  the  mistuning  problem,  researchers  have  developed 
reduced-order  models  (ROMs)  of  the  biaded  disk.  These  ROMs 
have  the  structural  fidelity  of  a  finite  element  model  of  the  full 
rotor,  while  incurring  computational  costs  that  are  comparable  to 
that  of  a  mass-spring  model,  [2-5].  In  numerical  simulations, 
most  published  ROMs  have  correlated  extremely  well  with  nu¬ 
merical  benchmarks.  However,  some  models  have  at  times  had 
difficulty  correlating  with  experimental  data,  [6].  These  results 
suggest  that  the  source  of  the  error  may  lie  in  our  inability  to 
determine  the  correct  input  parameters  to  the  ROMs. 

The  standard  method  of  measuring  mistuning  in  rotors  with 
attachable  blades  is  to  mount  each  blade  in  a  broach  block  and 
measure  its  natural  frequency.  The  difference  of  each  blade’s  natu¬ 
ral  frequency  from  the  mean  value  is  then  taken  as  a  measure  of 
the  mistuning.  However,  the  mistuning  measured  through  this 
method  may  be  significantly  different  from  the  mistuning  present 
once  the  blades  are  mounted  on  the  disk.  This  variation  in  mis¬ 
tuning  can  arise  because  each  blade’s  frequency  is  dependent  on 
the  contact  conditions  at  the  attachment.  Not  oiily  may  the  blade- 
broach  contact  differ  from  the  blade-disk  contact,  but  the  contact 
conditions  can  also  vary  from  slot-to-slot  around  the  wheel. 
Therefore,  in  order  to  accurately  measure  mistuning,  we  must  de¬ 
velop  methods  that  can  make  measurements  of  the  blade-disk  as¬ 
sembly  as  a  whole. 

Such  holistic  measurement  techniques  are  particularly  impor¬ 
tant  for  integrally  biaded  rotors,  since  their  blades  caimot  be  re¬ 
moved  for  individual  measurement.  In  this  paper,  we  present  a 
new  method  of  identifying  mistuning  in  biaded  disks  that  is  based 
on  the  vibration  characteristics  of  the  whole  system.  The  key  con- 
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cept  is  that  system  modes  are  highly  sensitive  to  small  variations 
in  mistiming.  Consequently,  the  modes  themselves  provide  a  sen¬ 
sitive  and  accurate  basis  for  identifying  mistuning  in  the  system. 

Oiir  method  of  system  identification  is  based  on  a  recently  de¬ 
veloped  reduced  order  model  called  the  fundamental  mistuning 
model  (FMM),  [5],  that  accurately  represents  the  vibratory  re¬ 
sponse  of  an  isolated  family  of  modes.  FMM  is  a  highly  reduced- 
order  model  that  can  completely  describe  a  mistuned  rotor  using 
only  its  tuned  system  frequencies  and  the  frequency  mistuning  of 
each  blade/disk  sector.  As  a  result,  when  the  FMM  based  identi¬ 
fication  method  (FMM  ID)  is  applicable,  it  is  very  easy  to  use  and 
requires  very  little  analytical  information  about  the  system,  e.g., 
no  finite  element  mass  or  stiffness  matrices.  We  have  developed 
two  forms  of  FMM  ED:  a  basic  version  of  FMM  ID  that  requires 
some  information  about  the  system  properties,  and  a  somewhat 
more  advanced  version  that  is  completely  ex^rimentally  based. 
The  basic  FMM  ID  requires  the  frequencies  of  the  tuned  system^ 
as  input.  Then,  given  measurements  of  a  limited  number  of  mis¬ 
tuned  modes  and  frequencies,^  FMM  ID  solves  for  the  mistimed 
frequency  of  each  sector.  The  advanced  form  of  FMM  ED  uses 
measurements  of  some  mistuned  modes  and  frequencies  to  deter¬ 
mine  ail  of  the  parameters  in  FMM,  i.e.,  the  frequencies  that  the 
system  would  have  if  it  were  tuned  as  well  as  the  mistuned  fre¬ 
quency  of  each  sector.  Thus,  the  tuned  system  frequencies  deter¬ 
mined  from  the  second  method  can  also  be  used  to  validate  finite 
element  models  of  the  nominal  system. 

Judge  and  Pierre  have  developed  an  approach  for  determining 
mistuning  in  IBRs,  [7],  that  also  uses  mistuned  modes  and  fre¬ 
quencies  to  infer  the  rotor’s  mistuning.  However,  it  is  based  on  a 
more  complicated  reduced-order  model  that  results  in  a  more 
complex  methodology,  and  requires  significantly  more  analyti¬ 
cally  generated  input  data  to  implement.  However,  it  is  not  inher¬ 
ently  limited  to  an  isolated  family  of  modes  as  is  FMM  ID.  Mi- 
gnolet  and  Rivas-Guerra  have  also  studied  mistuning 
identification,  [8,9].  Their  focus,  however,  was  on  discerning  the 
difference  between  mass  and  stifrhess  mistuning  in  an  isolated 
blade. 

This  paper  is  Part  I  of  a  two  part  investigation  and  develops  the 

^The  frequencies  of  the  tuned  system  are  typically  calculated  using  a  Mte  ele¬ 
ment  analysis  of  a  single  blade/disk  sector  with  cycUc  symmetric  boundary  condi¬ 
tions  applied  to  the  disk. 

^Thc  modes  required  in  FMM  ID  are  the  circumferential  modes  that  correspond  to 
the  tip  displacement  of  each  blade  around  the  wheel. 
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theory  of  the  FMM  ID  methods,  then  confirms  their  applicability 
with  numerical  test  cases.  The  second  paper,  Part  n,  examines  the 
application  of  the  methods  to  actual  hardware. 

This  paper  is  organized  as  follows.  Section  2  presents  the  basic 
FMM  ID  theory,  and  provides  numerical  test  cases  of  the  method. 
In  Section  3,  we  present  the  completely  experimental  version  of 
FMM  ID  that  also  identifies  the  tuned  system  frequencies.  Section 
4  discusses  an  extrapolation  method  that  allows  mistuning 
measured  at  rest  to  be  used  to  predict  the  response  of  a  rotor  at 
speed.  Lastly,  the  key  attributes  of  the  method  are  summarized  in 
Conclusions. 


2  Basic  FMM  ID 

This  section  presents  the  basic  FMM  ID  method.  The  basic 
method  uses  tuned  system  frequencies  along  with  measurements 
of  the  mistimed  rotor’s  system  modes  and  frequencies  to  infer 
mistuning, 

2.1  Theory.  The  FMM  ID  method  is  derived  from  the  fun¬ 
damental  mistuning  model  (FMM).  The  first  part  of  this  section 
describes  the  FMM  method.  Then  we  invert  the  equations  to  ob¬ 
tain  a  formulation  that  can  be  used  for  system  identification! 

2. LI  FMM  Method.  FMM  is  a  highly  simplified  reduced- 
order  model  that  can  accurately  predict  the  vibratoiy  response  of 
redistic  bladed  disks  in  an  isolated  family  of  modes,  [5].  FMM 
requires  only  two  sets  of  input  parameters  to  calculate  the  modes 
and  natural  frequencies  of  the  mistimed  system:  the  tuned  system 
frequencies  of  one  isolated  family  of  modes,  and  the  frequency 
deviation  of  each  blade-disk  sector.  In  Appendix  A,  we  derived  a 
more  general  form  of  the  FMM  modal  equation  than  was  given  in 
[5]  that  is  applicable  to  systems  with  more  flexible  disks.  The 
more  general  modal  equation  is 

.  (1) 

The  eigenvector  of  this  equation,  ySy,  contains  weighting  factors 
that  describe  the  yth  mistuned  mode  as  a  sum  of  tuned  modes,  i.e., 

N-\ 

(2) 

m=0 

where  is  the  mth  tuned  mode  of  the  family  of  interest.  The 
corresponding  eigenvalue,  (o],  is  the  ythmode’s  natural  frequency 
squared. 

Tl^  matrix  of  the  eigenvalue  problem  contains  two  terms,  11“ 
and  D.  XV  is  a  diagonal  matrix  of  the  tuned  system  frequencies, 
ordered  by  ascending  inter-blade  phase  angle  of  their  correspond¬ 
ing  mode.  The  notation  11“^  is  shorthand  for  0“^I1“,  which  results 
in  a  diagonal  matrix  of  the  tuned  system  frequencies  squared.  The 
matrix  11  contains  the  discrete  Fourier  transforms  (DFT)  of  the 
sector  frequency  deviations.  D  has  the  form 


«0 

"1  • 

®iV-I 

^0  * 

*  ^N-2 

. 

a>2  ' ' 

Wq  _ 

where  cop  is  the  pth  DFT  of  the  sector  frequency  deviations.  Note 
that  fl  is  a  circulant  matrix,  in  which  each  column  is  equal  to  the 
previous  column  rotated  down  a  row.  Therefore,  for  an  N  bladed 
disk,  it  has  only  N  distinct  values. 

One  of  the  key  changes  in  the  generalized  FMM  is  that  it  uses 
a  new  quantity  called  a  “sector  frequency  deviation”  as  a  measure 
of  the  mistuning  for  each  blade-disk  sector.  In  the  original  FMM 
formulation,  mistuning  was  measured  by  blade  frequency  devia¬ 
tions.  The  advantage  of  the  new  mistiming  measure  is  that  it  not 
only  accounts  for  mistuning  in  the  blade,  but  also  captures  mis- 
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tuning  in  the  disk  as  well  as  variations  in  the  way  the  blades  are 
attached  to  the  disk.  The  definition  of  a  sector  frequency  deviation 
is  provided  in  Appendix  A. 

FMM  treats  the  rotor’s  mistuning  as  a  known  quantity  that  it 
uses  to  determine  the  system’s  mistuned  modes  and  frequencies. 
However,  if  we  were  to  treat  the  mistuned  modes  and  frequencies 
as  known,  we  could  solve  the  inverse  problem  to  determine  the 
rotor’s  mistuning.  This  is  the  basis  of  FMM  ID. 

2.1.2  Inversion  of  FMM  Equation.  This  section  manipulates 
the  FMM  equation  of  motion  to  solve  for  the  mistuning  in  the 
mtor.  Consider  Eq.  (1).  All  quantities  are  treated  as  known  except 
fl,  which  describes  the  system’s  mistuning.  Subtracting  the 
term  from  both  sides  of  (1)  and  regrouping  terms  yields 

(4) 

The  bracketed  quantity  on  the  left-hand  side  of  (4)  contains  a 
known  vector,  which  will  be  denoted  as  yj , 

(5) 

Thus,  yj  simply  contains  the  modal  weighting  factors,  ySy ,  scaled 
on  an  element-by-element  basis  by  their  corresponding  natural 
frequencies.  Substituting  yj  into  (4)  yields 

2D“[fi7y]  =  .  (6) 

Consider  the  bracketed  term  of  this  expression.  After  some  al¬ 
gebra,  it  can  be  shown  that  this  product  may  be  rewritten  in  the 
form 

nyy-FyS  (7) 

where  the  vector  «  equals  The  matrix  Tj 

is  composed  from  the  elements  in  yj  and  has  the  form 
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where  yj^  denotes  the  nth  element  of  the  vector  yj ;  the  yj  ele¬ 
ments  are  numbered  from  0  to  iV-l. 

Substituting  (7)  into  (6)  produces  an  expression  in  which  the 
matrix  of  mistuning  parameters,  D,  has  been  replaced  by  a  vector 
of  mistuning  parameters,  6) 

(9) 

Observe  that  pre-multiplying  (9)  by  (2JIT^)"*  would  solve  this 
expression  for  the  DFT  of  the  rotor’s  mistuning.  Furthermore,  the 
vector  5}  can  then  be  related  to  the  physical  sector  mistuning 
through  an  inverse  discrete  Fourier  transfonn.  However,  (9)  only 
contains  data  from  one  measured  mode  and  frequency.  Therefore, 
error  in  the  mode’s  measurement  may  result  in  significant  error  in 
the  predicted  mistuning. 

To  minimize  the  effects  of  measurement  error,  we  v/ill  incorpo¬ 
rate  multiple  mode  measurements  into  our  solution  for  the  mis¬ 
tuning.  We  construct  (9)  for  each  of  the  M  measured  modes,  and 
combine  them  into  the  single  matrix  expression, 
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For  brevity,  we  rewrite  (10)  as 

La5  =  r  (11) 
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where  L  is  the  matrix  on  the  left-hand  side  of  the  expression,  and 
r  is  the  vector  oii  the  right-hand  side.  The  ”  is  used  to  indicate 
that  these  quantities  are  composed  by  vertically  stacking  a  set  of 
submatrices  or  vectors.  . 

Note  that  expression  (11)  is  an  over  determined  set  of  equa¬ 
tions.  Therefore,  we  can  no  longer  solve  for  by  direct  inverse. 
Howeyer,  we  can  obtain  a  least  squares  fit  to  the  mistuning,  i.e., 

B^Lsq{L,r},  (12) 

Equation  (12)  produces  the  vector  c5  that  best  fits  all  the  measured 
data.  Therefore,  the  error  in  each  measurement  is  compensated  for 
by  the  balance  of  the  data.  The  vector  w  can  ttien  be  related  to  the 
physical  sector  mistuning  through  the  inverse  transform, 

N-i 

(13) 

where  is  the  sector  frequency  deviation  of  the  ith  sector. 
The  following  section  describes  how  Eqs.  (12)  and  (13)  can  be 
applied  to  determine  a  rotor’s  mistuning. 

2JJ  Experimental  Application,  In  order  to  solve  Eqs.  (12) 

and  (13)  for  the  sector  mistuning,  we  must  first  construct  L  and  r 
from  the  tuned  system  frequencies  and  the  mistuned  modes  and 
frequencies.  The  tuned  system  frequencies  can  be  calculated 
through  finite  element  analysis  of  a  tuned,  cyclic  symmetric, 
single  blade/disk  sector  model.  However,  the  mistuned  modes  and 
frequencies  must  be  obtained  experimentally. 

The  modes  used  by  FMM  ID  are  circumferential  modes,  corre¬ 
sponding  to  the  tip  displacement  of  each  blade  on  the  rotor.  Since . 
EMM  ID  is  designed  for  isolated  families  of  modes,  it  is  sufficient 
to  measure  the  displacement  of  only  one  point  per  blade.  In  prac¬ 
tice,  modes  ^d  frequencies  are  obtained  by  first  measuring  a 
complete  set  of  frequency  response  functions  (FRFs).  Then,  the 
modes  and  frequencies  are  extracted  from  the  FRFs  using  modal 
curve  fitting  software. 

The  mistuned  frequencies  obtained  from  the  measurements  ap¬ 
pear  explicitly  in  the  FMM  ID  equations  as  ct>j  .*  However,  the 
mistuned  modes  enter  into  the  equations  indirectly  through  the 
modal  weighting  factors  pj .  As  described  by  Feiner  and  GriMn 
[5],  each  vector  is  obtained  by  taking  the  inverse  discrete  Four 
tier  transform  of  the  corresponding  single  point-per-blade  mode, 
i.e., 

N-l 

(14) 

in=0  . 

These  quantities  may  then  be  used  with  the  tuned  system 

frequencies  to  construct  L  and  f  as  outlined  in  earlier  portions  of 
this  section.  Finally,  (12)  and  (13)  may  be  solved, for  the  sector 
inistuning. 

This  process  is  denionstrated  through  the  two  examples  in  the 
following  section. 

2,2  Numerical  Examples.  This  section  presents  two  nu¬ 
merical  examples  of  the  basic  FMM  ID  metiiod.  In  the  first  ex¬ 
ample,  we  consider  an  integrally  bladed  compressor  whose  blades 
are  geometrically  mistuned.  The  sector  frequency  deviations  iden¬ 
tified  by  FMM  ID  are  verified  by  coriaparing  them  with  values 
directly  determined  by  finite  element  analyses.  The  second  test 
case  highlights  FMM  ID’s  ability  to  detect  mistuning  caused  by 
variations  at  the  blade-disk  interface. 

2.2,1  Geometric  Blade  Mistiming.  Consider  the  finite  ele¬ 
ment  model  of  the  twenty  blade  compressor  shown  in  Fig,  1. 
Although  the  airfoils  on  this  model  are  simply  flat  plates,  the  rotor 
design  reflects  the  key  dynamic  behaviors  of  a  modem,  integrally 
bladed  compressor.  We  mistuned  the  rotor  through  a  combination 
of  geometric  and  material  property  changes.  Approximately  one- 
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Fig.  1  Rnite  element  model  used  to  represent  compressor 


third  of  the  blad^  were  mistuned  through  length  variations,  one- 
third  through  thickness  variations,  and  one-third  through  elastic 
modulus  variations.  The  magnitudes  of  the  variations  were  chosen 
so  that  each  form  of  mistuning  would  contribute  equally  to  a  1.5% 
standard  deviati^on  in  the  sector  frequencies. 

We  first  performed  a  finite  element  analysis  of  the  tuned  rotor, 
and  generated  its  nodal  diameter  map.  Fig,  2.  Observe  that  the 
lowest  frequency  family  of  first  bending  modes  is  isolated,  and  is 
therefore  a  good  candi^te  for  FMM  ID.  The  sector  mistuning  of 
this  rotor  was  then  determined  through  two  different  methods: 
finite  element  analyses  of  the  mistuned  sectors  using  the  commer¬ 
cially  available  ANSYS  finite  element  code,  and  FMM  ID. 

The  finite  element  calculations  serve  as  a  benchmark  to  assess 
the  accuracy  of  the  FMM  ID  method.  In  the  benchmark,  a  finite 
element  mc^el  was  made  for  .each  mistuned  blade.  lii  the  model 
the  blade  is  attached  to  a  single  disk  sector.  The  frequency  change 
in  the mistuned  blade/disk  sector  was  then  calculated  with  various 
cyclic  symmetric  boundary  conditions  applied  to  the  disk.  It  was 
found  tl^t  file  phase  angle  of  tiie  cyclic  syminetric  constraint  had 
little  effect  on  the  frequency  change  caused  by  blade  mistuning. 
The  values  quoted  in  this  paper  are  for  a  disk  phase  constraint  of 
90  deg,  i.e.,  for  the  five  nodal  diameter  mode. 


Fig.  2  Natural  frequencies  of  the  compressor  with  no 
mistuning 
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A  finite  element  model  of  the  full,  mistuned  bladed  disk  was 

also  constructed  and  used  to  compute  its  mistuned  modes  and  variation  of  ±4  deg.  Otherwise  the  blades  are  identical.  The 

natural  firequencies.  The  modes  and  fiequencies  were  used  as  in-  ^odes  of  the  system  were  then  calculated  using  the  ANSYS  finite 

put  data  for  FMM  ID»  In  an  actual  experiment,  the  mistuned  element  code,  figure  5  shows  a  r^resentativc  mode.  Observe 

modes  and  firequencies  would  be  obtained  through  a  modal  fit  of  localized,  indicating  that  varying  the  stagger 

the  rotor’s  fi^uency  response  functibns.  Typically,  die  measure-  angles  does  indeed  mistune  the  system. 

ments  will  not  det^  modes  that  have  a  node  point  at  the  excita-  mistuned  mbdes  and  fir^uencies  calculated 

tion  source.  To  reflect  this  phenomenon  in  our  numerical  test  case,  ANSYS  to  perform  an  FMM  ID  analysis  of  the  mistuning.  The 
we  eliminated  all  mistuned  modes  that  had  a  sinaH  response  at  resulting  sector  frequency  deviations  are  plotted  as  the  solid  line 

blade  one.  This  left  16  modes  and  natural  firequencies  to  apply  to  ^  The  squares  correspond  to  the  stagger  angle  variations 

FMM  ID.  applied  to  each  blade.  The  vertical  axes  have  been  scaled  so  that 

The  mistuned  modes  and  frequencies  were  combined  with  the  maximum  frequency  and  angle  variation  data  points  (blade  14) 
tuned  system  frequencies  of  the  fundamental  mode  family  to  con-  ^  coincident.  This  was  done  to  highlight  the  fact  fiiat  the  stagger 

struct  the  basic  FMM  ID  equations,  (10)^  These  equations  were  variations  arc  proportional  to  the  sector  frequency  devia- 

solved  using  a  least-squ^s  fit  The  solution  was  then  converted  detected  by  IMM  ID.  Thus,  not  only  can  FMM  ID  accu- 

to  the  physical  sector  frequency  deviations  through  the,  inverse  retely  detect  mistuning  in  the  blades,  as  illustrate  in  the  previous 
transform  given  in  (13),  example^  but  it  can  also  accurately  detect  other  foimis  of  mistun- 

Fignre  3  shows  the  comparison  between  the  sector  mistuning  such  as  variation  in  the  blade  stagger  angle, 
calculated  directly  by  finite  element  simulations  of  each  mistuned 

blade/sector  and  the  mistuning  identified  by  FMM  JD.  The  two  3  A  Completely  Experimental  Method  ot  Identifica* 
results  are  in  good  agreement.  tion 
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technique  requires  a  set  of 'simple  vibration  measurements  and  the 
'  natural  firequencies  of  the  tuned  system.  However,  at  times  neither 
the  tuned  system  frequencies  nor  a  finite  element  model  from 
which  to  obtain  them  are  available  to  researchers  interested  in 
determining  an  IBR’s  mistuning.  Furthermore,  even  if  a  finite  el¬ 
ement  model  is  available,  there  is  often  concern  as  to  how  accu¬ 
rately  the  model  represents  the  actual  rotor.  Iherefore,  we  have 
developed  an  alternative  FMM  ID  method  that  does  not  require 
any  analytical  data.  The  approach  requires  only  a  limited  number 
of  mistuned  modes  and  frequency  measurements  to  determine  a 
bladed  disk’s  mistuning.  Furthermore,  the  method  also  identifies 
the  bladed  disk’s  tuned  system  frequencies.  Thus,  it  not  only 
serves  as  a  method  of  identifying  mistuning  the  system,  but  can 
also  provide  a  method  of  corroborating  the  finite  element  model  of 
.  the  tuned  system. 

3.1  Theory.  This  version  of  FMM  ID  is  derived  from  the 
basic  FMM  ID  equations.  Recall  that  an  important  step  in  the 
development  of  the  basic  FMM  ID  theory  was  to  transform  the 
mistuning  matrix  £1  into  a  vector  form.  Once  we  expressed  the 
mistuning  as  a  vector,  it  could  then  be  calculated  using  standard 
methods  from  linear  algebra.  A  similar  approach  is  used  in  the 
current  development  to  solve  for  the.  tuned  system  frequencies. 
However,  the  resulting  equations  are  nonlinear,  and  require  a  more 
sophisticated  solution  approach. 

5.7./  Development  of  Nonlinear  Equations.  Consider  the  ba¬ 
sic  FMM  ID  equation  given  in  (9).  Moving  the  £1*^  term  to  the 
left-hand  side,  the  expression  becomes 

(15) 


We  assume  that  from  measurement  of  the  mistuned  modes  and 
frequencies,  and  (Oj  are  known.  AH  other  quantities  are  un? 
known.  Note  that  although  Fj  is  not  known,  the  matrix  contains 
elements  from  pj .  Therefore,  we  do  have  some  knowledge  of  the 
matrix. 

After  some  algebra,  one  can  show  that  the  term  may  be 
re-expressed  as 

(16) 


where  a  vector  of  the  tuned  frequencies  squared,  and  Bj  is  a 
matrix  composed  from  the  elements  of  Pj .  If  we  define  ^7  to  be 
the  maximum  number  of  nodal  diameters  on  the  rotor,  i.e.,  77=N/2 
if  N  is  even  or  (N-I)/2  if  N  is  odd,  then  X*  is  given  by 
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For  N  even,  the  matrix  has  the  form 
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Equation  (19)  contains  information  from  only  one  of  the  M 
measured  modes  and  frequencies.  However,  (19)  can  be  con¬ 
structed  for  each  nieasured  mode,  and  combined  into  the  single 
matrix  expression 


(20) 


Thus  we  have  formed  a  single  expression  that  incorporates  all  of 
the  measured  data.  For  brevity,  (20)  is  rewritten  ^ 
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where  B  is  the  stacked  matrix  of  the  term  (fFF)  is  the 

stacked  matrix  of  ^FFj ,  and  f '  is  the  right-hand  side  of  (20). 

To  complete  our  formulation,  we  must  introduce  ah  additional 
constraint  equation.  It  is  required  because  the  Eqs.  (21)  are  under- 
determin^  To  understand  the  cause  of  diis  indeterminacy,  con¬ 
sider  a  rotor  in  which  each  sector  is  mistuned  the  same  amount 
Due  to.  the  syinmetry  of  the  mistuning,  the  rotor’s  mode  shapes 
will  still  look  tuned,  but  its  frequencies.  wiQ  be  shifted.  If  one  has 
no  prior  knowledge  of  the  tuned  system  frequencies,  there  is  no 
way  to  determine  that  the  rotor  has  in  fact  been  mistuned.  The 
same  difficulty  arises  in  solving  (21)  since  there  is  no  way  to 
distinguish  between  a  mean  shift  in  the  mistuning  and  a  corre¬ 
sponding  shift  in  the  tuned  system  frequencies.  To  eliminate  this 
ambiguity,  we  will  define  mistuning  so  that  it  has  a  mean  value  of 
zero. 

'  Mathematically,  a  zero  mean  in  the  mistuning  translates  to  pre¬ 
scribing  the  first  element  of  55  to  be  zero.  With  the  addition  of  this 
constraint,  (21)  takes  the  form 


(22) 


where  c  is  a  row  vector  whose  first  element  is  1  and  whose  re¬ 
maining  elements  are  zero. 
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5.7.2_  Iterative  Solution  Method,  Consider  Eq.  (22).  If  the 
term  (FFF)  were  known,  then  we  could  obtain  a  least-squares 
solution  for  tiie  tuned  eigenvalues^*  and  the  DFT  of  the  sector 
mistuning  c5.  However,  since  (£1*1)  is  based  in  part  on  the  un¬ 
known  quantities  X*,  tiie  equations  are  nonlinear.  Therefore,  we 
must  use  an  alternative  solution  method.  This  section  describes 
how  these  equations  may  be  solved  using  an  iterative  approach. 

In  iterative  form,  the  least  squares  solution  to  (22)  cmi  be  writ¬ 
ten  as 


x* 

=Li^| 
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2(«T)(,_,)' 

.a 

(i)  1 
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c  . 

(23) 


where  the  subscripts  indicate  the  iteratioj^number.  For  each  itera¬ 
tion,  we  then  construct  a  new  matrix  (£l*r)  based  on  the  previous 
iteration’s  solution  for  X*.  This  process  is  repeated  until  we  obtain 
a  converged  solution.  With  a  good  initial  guess,  this  method  typi¬ 
cally  conveiges  within  a  few  iterations.  Appendix  B  describes  an 
effective  approach  for  obtaining  a  good  initial  guess. 


A  similar  expression  can  be  derived  for  N  odd. 

Substituting  (16)  into  (15)  and  regrouping  the  left-hand  side 
results  in  a  matrix  equation  for  the  tuned  frequencies  squared  and 
the  sector  mistuning, 
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3.2  Numerical  Test  Case.  This  section  presents  a  numerical 
example  of  the  FMM  ID  method  that  identifies  the  tuned  system 
frequencies  as  well  as  the  mistuning.  This  example  uses  the  geo¬ 
metrically  mistuned  compressor  model  presented  in  Section  2. 
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4  R^ponse  Prediction  at  Speed 

The  FMM  ID  methods  presented  in  Sections  2  and  3  detennine 
the  mistuhing  in  bladed  disks  While  they  are  stationary.  However, 
once  the  rotor  is  spinning,  centrifugal  forces  can  alter  the  effective 
mistuning  in  the  rotor.  Section  4.1  presents  a  method  for  approxi¬ 
mately  calculating  the  effect  of  rotational  speed  on  mistuning. 
Then,  in  Section  4.2,  we  present  a  numerical  example  that  applies 
•  this  technique  and  then  uses  FMM  to  predict  its  forced  response  at 
speed.  The  accuracy  of  the  method  is  assessed  by  comparing  the 
results  with  a  numerical  benchmarL 

4 J  Mistuning  Extrapolation*  Centrifugal  effects  cause  the 
sector  frequency  deviations  to  change  when  the  disk  is  rotating. 
This  section  gives  a  simple  method  for  approximating  the  effect  of 
rotational  speed  on  mistuning.  Details  of  the  derivation  arc  pro¬ 
vided  by  Feiner  [10],  Here  we  summarize  the  result 

To  investigate  centrifugal  stiffening  in  [10]  the  blade  is  mod¬ 
eled  as  a  pendulum,  constrained  to  a  rigid  disk  by  a  torsion  spring. 
While  the  disk  is  at  rest,  the  blade’s  frequency  is  perturbed  (mis- 
tuned).  Then,  the  model  is  used  to  determine  the  corresponding 
change  in  the  blade’s  frequency  when  the  disk  is  rotating.  To  first 
order  one  finds  that 


The  tuned  system  frequencies  and  sector  mistuning  identified  by 
FMM  ID  are  then  compared  with  finite  element  results. 

We  calculated  the  modes  and  natural  frequencies  of  the  mis- 
tuned  bladed  disk  using  a  finite  element  model  of  the  mistuned 
system.  We  then  converted  the  physical  modes  to  vectors  of  modal 
weighting  fectors,  ,  through  (14).  The  weighting  factors  were 
used  to  form  the  elements  of  Eq.  (48)  in  Appendix  B  which  was 
solved  to  obtain  an  initial  estimate  of  the  umed  system  frequen¬ 
cies.  This  was  used  as  an  initial  guess  to  iteratively  solve  Eq.  (23). 
The  solution  vector  contains  two  parts:  a  vector  of  the  tuned  sys¬ 
tem  frequencies  squared,  and  a  vector  of  the  DFT  of  the  sector 
frequency  deviations.  The  sector  mistuning  was  converted  to  the 
physical  domain  using  the  inverse  transform  (13). 

The  resulting  sector  frequency  deviations  are  compared  with 
the  benchmark  finite  element  values  in  Fig.  7  using  the  same 
procedure  as  in  Section  2.2.  A  comparison  of  the  tuned  frequen¬ 
cies  identified  by  FMM  ID  and  those  computed  directly  with  the 
finite  element  model  is  shown  in  Fig.  8.  In  each  case  the  agree¬ 
ment  is  good. 


Rg.  8  Comparison  of  the  tuned  system  frequencies  from  FMM 
ID  and  FEM 


Acti(5)<'>-A<w(0)^"> 


(24) 


where  S  is  the  rotation  speed  of  the  disk,  Aco^S)^^^  is;  the  sector 
mistuning  ratio  of  the  sth  sector  and  ci>^^(5)  is  the  frequency  of 
a  representative  timed  system  mode.  For  example,  in  the  case 
r^rted  in  4.2,  6>^p(5)  is  the  frequency  of  the  system  mode 
with  a  90  deg  interblade  phase  angle.  A  key  result  from  (24)  is  that 
mistuning  is  larger  at  rest  than  at  speed,  and  that  this  reduction  in 
mistuning  can  be  estimated  by  calculating  how  the  tuned  system 
frequencies  change  with  speed.  It  has  been  confirmed  that  (24)  - 
works  reasonably  well  for  realistic  geometric  mistuning  in  real 
compressor  blades,  [10].  This  will  be  also  be  demonstrated  by  the 
next  example. 

4.2  Numerical  Test  Case.  This  section  uses  a  numerical 
test  case  that  shows  how  FMM  ID,  Eq.  (24),  and  the  FMM  forced 
response  code  can  be  combined  to  predict  the  response  of  a  bladed 
disk  under  rotating  conditions.  Consider  the  geometrically  mis- 
tuned  rotor  studied  in  Section  2.2.1,  and  illustrated  in  Fig.  1.  This 
compressor  has  a  6th  engine  order  crossing  ,  with  the  first  bending 
modes  at  a  rotational  speed  of  20,000  rpm.  However,  to  create  a 
more  severe  test  case,  we  will  proceed  as  if  the  crossing  occurred 
at  40,000  rpm. 

In  order  to  use  FMM  to  predict  the  rotor’s  forced  response  at 
this  speed,  we  must  provide  the  FMM  prediction  code,  [5],  with 
the  bladed  disk’s  tuned  system  frequencies  and  the  sector  fre¬ 
quency  deviations  that  are  present  at  40,000  rpm,  La  Section  2.2, 1 
we  determined  these  two  sets  of  parameters  at  rest  using  ANSYS 
and  Basic  FMM  ID,  respectively.  However,,  since  both  of  these 
properties  change  witii  rotation  speed,  they  must  first  be  adjusted 
to  reflect  their  values  at  40,000  ipm. 

To  adjust  die  tuned  system  frequencies,  we  recalculated  tfaeni  in 
ANSYS  using  the  centrifugal  load  option  to  simulate  rotational 
effects.  The  centrifugal  stiffening  caused  the  tuned  system  fre¬ 
quencies  to  increase  by  about  30%.  Then  we  used  the  change  in 
the  five  nodal  diameter,  tuned  system  frequency  and  (24)  to  ana¬ 
lytically  extrapolate  the  sector  frequency  deviations  to  40,000 
rpm.  In  this  case,  the  centrifugal  loading  reduces  the  mistuning 
ratios  by  about  40%. 

Finally,  the  adjusted  parameters  were  used  with  the  FMH 
forced  response  code  to  calculate  the  rotor’s  response  to  a  6E 
excitation  using  the  method  described  in  [5].  As  a  benchmark,  the 
forced  response  was  also  calculated  directly  in  ANSYS  using  a 
foil  360  deg  mistuned  finite  element  model.  Trackirig  plots  of  the 
FMM  and  ANSYS  results  are  shown  in  Hg.  9,  For  clarity,  we 
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Excitation  Frequency  (Hz) 

Fig.  9  Frequency  response  of  blades  to  a  six  engine  order 
excitation 


have  shown  the  response  of  only  three  blades:  the  high  responding 
blade,  the  median  responding  blade,  and  the  low  responding 
blade.  Observe  that  each  blade’s  peak  amplitude  and  the  shape  of 
its  overall  response  as  predated  by  FMM  agree  well  with  the 
benchmark  results.  Thus,  by  combining  FMM  ID,  the  mistuning 
extrapolation  equation,  and  FMM,  we  were  able  to  identify  the 
mistuning  of  a  rotor  at  rest,  and  use  it  to  accurately  predict  the 
system’s  forced  response  imder  rotating  conditions. 

5  Conclusions 

A  new  method  of  identifying  mistuning  in  bladed  disks  is  de¬ 
veloped-  This  approach  is  called  FMM  ID  because  it  is  based  on 
the  fundamental  mistuning  model  (FMM).  FMM  and  the  methods 
presented  in  this  paper  are  applicable  to  isolated  families  of 
modes.  Often  the  frequencies  of  the  first  bending  and  first  torsion 
families  of  modes  satisfy  this  requirement  Identifying  mistuning 
in  these  inodes  is  important  not  only  for  predicting  forced  re¬ 
sponse,  but  also  for  predicting  flutter.  Since  mistuning  tends  to 
stabilize  flutter,  tiie  issue  is  how  to  relate  a  particular  flutter  test  to 
the  fleet  as  a  whole.  We  are  pursuing  the  application  of  FMM  ID 
to  the  flutter  problem  in  a  joint  research  effort  with  Honeywell 
research  engineers. 

FMM  ID  uses  measurements  of  the  system  mode  shapes  and 
natural  frequencies  to  infer  the  rotor’s  mistuning.  The  key  concept 
behind  FMM  ID  is  that  the  high  sensitivity  of  system  modes  to 
small  variations  in  mistuning,  causes  measurements  of  those 
modes  themselves  to  be  an  accurate  basis  for  mistuning  identifi¬ 
cation.  Since  FMM  ID  does  not  require  individual  blade  measure¬ 
ments,  it  is  particularly  suited  to  integrally  bladed  rotors.  We  have 
developed  two  forms  of  FMM  ID:  Basic  FMM  ID  and  a  version 
that  also  identifies  the  frequencies  of  the  tuned  system. 

Basic  FMM  ID  uses  tuned  system  frequencies  from  finite  ele¬ 
ment  analysis  and  measurements  of  the  mistuned  system  modes 
and  frequencies  to  determine  a  rotor’s;  mistuning.  The  mistuned 
mckie  slmpe  and  frequencies  can  be  measured  with  standard  modal 
analysis  techniques. 

Since  Basic  FMM  ID  is  derived  from  FMM,  it  requires  very 
little  analytical  information.  Specifically,  it  needs  only  the  tuned 
system  frequencies  of  the  mode  family  of  interest  Thus,  for  an  N 
bladed  disk,  we  only  need  approximately  N/2  pieces  of  analytical 
data.  These  frequencies  can  be  obtained  from  finite  element  analy¬ 
sis  of  a  single-sector  model  using  cyclic  symmetric  boundary  con¬ 
ditions.  The  alternate  form  of  FMM  ID  requires  no  analytical  data. 
It  relies  solely  on  experimental  measurements  of  the  mistuned 
modes  and  frequencies.  Thus,  the  second  form  of  FMM  ID  can  be 
u^d  to  identify  mistuning  even  if  a  finite  element  model  of  the 
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bladed  disk  is  not  available:  Furthermore,  this  approach  not  only 
identifies  the  mistuning  m  a  n>tor,  but  it  also  infers  its  tuned  sys¬ 
tem  frequencies.  Identifying  the  tuned  system  frequencies  may  be 
particularly  useful  for  assessing  the  v^dity  of  a  finite  element 
model  of  the  nominal  system. 

A  number  of  numeiicai  test  cases  are  analyzed  to  demonstrate 
the  applicability  of  the  methods.  One  of  these  involve  introducing 
mistuning  by  varying  the  stagger  angle  of  each  blade  in  what  w^ 
otherwise  a  perfectly  tuned  system.  FMM  ID  accurately  detects 
the  pattern  of  the  stagger  angle  mistuning.  Tlus  example  is  impor¬ 
tant  because  it  illustrates  the  fact  that  mistuning  in  the  bladed  disk 
can  be  caused  by  sources  that  cannot  be  measured  simply  in  terms 
of  blade  frequencies. 

FMM  ID  can  be  used  to  identify  the  mistuning  in  a  bladed  disk 
when  it  is  tested  in  the  laboratory.  A  method  has  been  demon¬ 
strated  for  ^proximating  how  centtifugal  loading  will  change  the 
mistuning  when  it  is  rotating  in  the  engine.  Other  factors  may  also 
be  present  in  the  engine  that  can  affect  the  mistun^  response. 
These  may  include:  temperature  effects,  gas  bending  stresses,  how 
the  disk  is  constrained  in  the  engine,  and  how  the  teeth  in  the 
attachment  change  their  contact  if  the  blades  are  conventionally 
attached  to  the  disk.  Except  for  the  constraints  on  the  disk,  these 
additional  effects  may  be  relatively  unimportant  in  integrally 
bladed  compressor  stages.  The  disk  constramts  can  be  taken  into 
account  by  perfomung  the  system  ID.  on  the  IBR  after  the  ftiU 
rotor  is  assembled.  Consequently,  it  seems  feasible  that  the  metii- 
odology  presented  in  this  paper  can  be  used  to  predict  the  vibra¬ 
tory  response  of  actual  compressor  stages:  U.S.  Air  Force  engi¬ 
neers  have  agreed  to  try  this  approach  in  forthcoming  tests  in  the 
Compressor  Research  Facility  (CRF)  at  Wright  Patterson  Air 
Force  Base  in  Dayton,  Ohio.  If  the  approach  proves  successful, 
then  the  plan  is  to  use  the  methodology  to  select  which  blades  will 
be  instrumented,  interpret  test  data,  and  relate  the  vibratory  re¬ 
sponse  measured  in  the  CRF  to  the  vibration  that  will  occur  in  the 
fleet  as  a  whole  . 
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Appendix  A  . 

A  More  GJeneral  Form  of  FMM.  This  appendix  presents  the 
derivation  of  a  more  general  form  of  the  modal  equation  for  the 
fundamental  mistuning  model  (FMM)  that  is  applicable  to  rotors 
mth  more  flexible  disks.  The  generalized  FMM  formulation  dif¬ 
fers  from  the  original  in  two  ways.  First,  it  no  longer  approxi¬ 
mates  tuned  system  frequencies  by  their  average  value.  This 
allows  for  a  much  larger  variation  among  the  tuned  frequencies. 
Second,  rather  than  using  the  blade-alone  mode  as  an  approxima¬ 
tion  of  the  various  nodal  diameter  sector  modes,  we  now  use  a 
representative  mode  of  a  single  blade-disk  sector.  Consequently, 
the  approach  now  includes  the  disk  portion  of  the  mode  shape, 
and  thus  allows  for  mote  strain  energy  in  the  disk. 

The  changes  in  the  formulation  also  modify  our  measure  of 
mistuning.  In  the  original  FMM  form,  we  measured  mistuning  as 
a  percent  deviation  in  the  blade-alone  frequency.  However,  mis¬ 
tuning  IS  now  measured  as  a  percent  deviation  the  frequency  of 
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sector;  The  advantage;!^  that  the  sector  frequency 
-deviations  not  only  capture  mistuning  in  the  hlade,  but  can  also 
capture  mistuning  in  the  disk  as  well  as  variations  in  the  ways  the 
blades  are  attached  to  the  disk. 

Section  1  of  this  appendix  describes  how  the  SNM  approach, 
11,  is  used  to  reduce  the  order  of  the  mistuned  free-response  equa¬ 
tions^  and  formulates  the  problem  in  terms  of  reduced-order  sector 
matrices.  Section  2  relates  the  sector  matrices  to  mistuned  sector 
frequencies.  Section  3  simplifies  the  resulting  mathematical  ex¬ 
pressions. 

1  Reduction  of  Order .  Consider  a  mistuned,  bladed  disk  in 
the  absence  of  an  excitation.  The  order  of  its  equation  of  motion  is 
reduced  through  a  subset  of  nominal  modes  approach.  The  result¬ 
ing  reduced-order  equation  can  be  written  as,  11: 

[(Xl*^+AK)““a>y(I+AM)]^y=0  (25) 

ft  ^  is  a  ^agonal  matrk  of  the  tuned  system  eigenvalues,^  and  I  is 
the  identity  matrix.  AK  and  AM  are  the  variations  in  the  modal 
stiffness  and  modal  m^s  matrices  caused  by  stiffness  and  mass 
mistuning.  The  vector  pj  contains  weighting  factors  that  describe 
the  Jih  mistuned  mode  as  a  limited  sum  of  tuned  modes,  i.e., 

(26) 

where  is  a  matrix  whose  columns  are  a  limited  number  of  the 
tuned  system  modes. 

Note  that  to  first  order,  (I+A]^->«(I-  AM).  Hius,  by  pre- 

muitiplying  (25)  by  (I+AIW[)“*  and  keeping  only  first-order 
terms,  the  expression  becomes 

(27) 

where 

A=  AK-  AMft*2.  (28) 

Next,  we  will  relate  the  matrix  A  to  the  frequency  deviations  of 
the  mistuned  sectors. 

2  Relating  Mistuning  to  Sector  Frequency  Deviations.  Re- 
la^g  A  to  frequency  deviations  is  a  three-step  process.  Hrst,  the 
mistumng  matrix  is  expressed  in  terms  of  the  system  noiode  shapes 
of  an  individual  sector.  Then,  the  system  sector  modes  are  related 
to  the  corresponding  mode  of  a  single,  isolated  sector.  Finally,  the 
resulting  sector-mode  terms  in  A  are  expressed  in  terms  of  the 
frequency  deviations  of  the  sectors. 

2.i  Relatbig  mistumng  to  system  sector  modes.  Consider  the 
mistuning  naatrix,  A,  in  (28).  This  matrix  can  be  expressed  as  a 
sum  of  the  contributions  from  each  mistuned  sector. 

N-\ 

A=2)  AW  (29) 

where  the  superscript  denotes  that  the  mistuning  corresponds  to 
the  jth  sector.  The  expression  for  a  single  element  of  A<^^  is 

6);2am<*))^;«  oo) 

where  AK^''^  and  AM^'*^  are  the  physical  stiffness  and  mass  per¬ 
turbations  of  the  jtfa  sector.  The  modes  and  are  the 
portions  of  the  mth  and  nth  columns  of  which  describe  the  sth 
sector  s  motion.  The  term  is  the  nth  diagonal  element  of 
Equation  (30)  relates  the  mistuning  to  the  system  sector  modes. 

In  the  next  section,  fliese  modes  are  related  to  the  mode  of  a  single 
isolated  blade-disk  sector. 

An  eigenvalue  is  equal  to  the  square  of  the  natural  frequency  of  a  mode. 
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1.^:2  Relating  system  sector  mod^  to  ah  average  iecior  mode 
The  tuned  modes  in  (30)  are  expressed  in  complex  traveling  wave 
form.  Thus,  the  motion  of  the  .rth  sector  can  be  related  to  the 
motion  of  the  0th  sector  by  a  phase  shift  This  allows  us  to  restate 
(30)  as 

Because  the  tuned  modes  used  in  flie  SNM  formulation  are  an 
isolated  family  of  modes,  the  sector  modes  of  afl  nodal  dimeters 
look  nearly  identical.  Therefore,  we  can  approximate  the  various 
sector  modes  by  an  average  sector  mode.  Applying  die  avenge 
^tor  mode  approximation  for  the  system  sector  modes  in  f31) 
A^^  can  be  written  as 

\  0)^  / 

(32) 

where  is  the  average  tuned  system  sector  mode,  and  is  its 
natural  frequency.  In  practice,  can  be  taken  to  be  the  median 
modal  diameter  mode.  The  factor  (w>;)/(<w^^)  scales  the  aver¬ 
age  sector  mode  terms  so  that  they  have  the  approximately  the 
same  strain  aicrgy  as  the  sector  modes  they  replace. 

2.3  Introduction  of  sector  frequency  deviation.  This  version 
of  FMM  uses  the  deviation  in  a  sector  frequency  quantity  to  mea¬ 
sure  mistuning.  To  understand  this  concept,  consider  an  imaginary 
“test”  rotor.  In  the  test  rotor  every  sector  is  inistuned  in  the  same 
fchion,  so  as  to  match  the  mistuning  in  our  sector  of  interest. 
Since  our  test  rotor's  mistuning  is  cyclically  symmetric,  its  mode 
shapes  are  vitally  identical  to  those  of  the  tuned  system.  How¬ 
ever,  there  wOI  be  a  shift  in  the  tuned  system  frequencies.  For 
small  levels  of  mistuning,  tiie  frequency  shift  is  nearly  the  same  in 
all  of  the  tuned  system  modes  and  can  be  approximated  by  the 
fractional  change  hi  the  frequericy  of  the  median  nodal  diameter 
mode.  Thus,  the  fractional  shift  in  the  median  nodal  diameter’s 
frequency  is  taken  as  our  measure  of  mistuning  and  is  defined  as 
the  sector  frequency  deviation. 

The  bracketed  terms  of  (32)  are  related  to  these  frequency  de¬ 
viations  ill  the  foUowing  manner.  Consider  a  bladed  disk  that  is 
mistuned  in  a  cyclic  symmetric  fashion,  i.e.,  each  sector  under¬ 
goes  the  same  mistuning.  Its  free-response  equation  of  motion  is 
given  by  the  expression 

[(K*-l-AK)-w^(]Vr-bAM)]^„=0.  (33) 

T^e  fte  mode  to  be  the^mistuned  version  of  the  tuned 
median  nodal  diameter  mode,  is  the  full  system  mode 

counterpart  of  the  average  sector  mode  Since  mistuning  is 
S3mimetric,  the  tuned  and  mistuned  versions  of  the  mode  are 
nearly  identical.  Substituting  for  and  pre-multiplying  by 
yields 

+  (34) 

These  terms  may  be  rearranged  to  isolate  the  frequency  terms 

^^^(AK-c.>2AM)^J’=6)?-iu;^  (35)  . 

Since  the  mistuning  is  symmetric,  each  sector  contributes  equally 
to  (35).  Thus,  the  contribution  from  the  0th  sector  is 

^«(0)ff(AK-a,2AM)^”W=i(n,2-n,;^).  (36) 

By  factoring  the  frequency  terms  oil  the  right-hand  side  of  (36),  it 

can  be  shown  that 

is  the  case  for  an  isolated  family  of  modes  in  ^diich  the  strain  enetgy  is 
primarily  in  the  blades.  If  there  is  a  significant  amount  of  strain  energy  in  the  disk 
then  the  firequency  of  the  modes  change  significantly  as  a  function  of  nodal  diameter 
and  the  modes  are  not  isolated,  i.e.,  they  cover  suc^  a  broad  frequeiuy  range  that 
they  interact  with  other  families  of  modes. 
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where  Aa>^  is  the  fractional  change  in  ^’s  natural  frequency  due 
to  mistuning,  given  by  Note  that  by  defi¬ 

nition  Aft)^  is  a  sector  frequency  deviation,  ^nation  (37)  can  be 
substituted  for  the  bracketed  terms  of  (32),  resulting  in  an  expres¬ 
sion  that  relates  the  elements  of  the  sector  s  mistuning  matrix  to 
that  sector’s  jfreqUency  deviation 

(38) 

where  the  superscript  on  Ao)^  is  introduced  to  indicate  that  the 
frequency  deviation  corresponds  to  the  5th  sector.  These  sector 
contributions  may  be  summed  to  obtain  the  elements  of  the  mis¬ 
tuning  matrix 

A„„=2a.>;[^S  .  (39) 

iV  ,=o 


■X(oj— Wavg  ■  '  '  "■('I4) 

However,  the  value  of  <«)*yg  is  not  known  and  therefore  cannot 
be  directly  applied  to  (23).  Consequently,  we  vwH  formulate  a 
slightly  modified  form  of  (22),  which  incorporates  the  initial  guess 
defined  by  (44^  Consider  Eq.  (22).  If  wejake  the  tuned  frequen¬ 
cies  to  be  equal  to  ,  tiien  the  term  (0°r)  may  be  expressed  as 

(St)=«;„^  (45) 

where  f  is  the  matrix  formed  by  vertically  stacking  the  MTj 
■  matrices. 

The  matrix  T;  is  also  related  to  the  tuned  frequencies.  As  a 
result,  the  elements  of  each  matrix  r^-  simplify  to  the  form 
■  This  allows  us  to  rewrite  as 

r,=  a.;vgZ,-  (46) 

where  Zj  is  composed  of  the  elements  arranged  in  the  same 
pattern  as  the  yj^  elements  shown  in  (8).  Thus,  consolidating  all 
terms,  (45)  can  be  written  as 


3  The  Simplified  Form  of  the  Fundamental  Mistuning  Model 
Modal  Equation,  The  bracketed  term  in  (39)  is  the  discrete  Fou¬ 
rier  transform  (DFT)  of  the  sector  frequency  deviations.  If  we  use 
the  dummy  variable  p  to  replace  the  quantity  (ri— m)  in  (39),  then 
the  pth  DFT  of  the  sector  frequency  deviations  is  given  by 

«P=  TjS  (40) 

’’NfTa. 


(arr)=(o2^  (47) 

w’here  Z  is  the  stacked  form  of  the  Zj  matrices. 

Substituting  (47)  into  (22)  and  regrouping  terms  results  in  the 
expression 


where  55^  denotes  the  pth  DFT.  By  substituting  (40)  into  (39),  A 
may  be  expressed  in  the  simplified  matrix  form 

.  (41) 

where 

■  So  ^N-l 

a.  ®“  ■■■  ®r  .  <«) 

_  Si  ©2  ***  . 

is  a  matrix  which  contains  die  discrete__Fourier  transfonns  of 
the  sector  frequency  deviations.  Note  that  has  a  circulant  form, 
and  thus  contains  only  N  distinct  elements,  fl’  is  a  diagonal  ma¬ 
trix  of  the  tuned  system  frequencies. 

Substituting  (41)  into  (27)  produces  the  most  basic  form  of  the 
eigenvalue  problem  that  may  be  solved  to  determine  the  modes 
and  natural  frequencies  of  the  mistuned  system. 

(43) 

Appendix  B 

Estimating  the  Thned  System  Frequencies.  This  appendix 
presents  an  effective  method  of  obtaining  a  good  initial  guess  of 
the  tuned  system  frequencies  for  use  in  the  iterative  solution  pro¬ 
cess  described  in  Section  3.1.2.  The  approach  is  to  obtain  the 
initial  guess  by  solving  a  companion  problem. 

To  identify  a  good  initial  guess,  recall  that  FMM  ID  requires 
that  we  analyze  an  isolated  family  of  modes.  In  general,  the  fre¬ 
quencies  of  isolated  mode  families  tend  to  span  a  fairly  small 
range.  Therefore,  they  may  be  reasonably  well  approximated  by 
their  mean  value,  i.e.. 


Note  that  the  term  was  grouped  yrith  the  vector  6).  Thus,  all 
the  unknown  expressions  are  consolidated  into  the  single  vector, 
on  the  left-hand  side  of  (48).  These  quantities  can  be  solved 
through  a  least  squares  fit  of  the  equations.  This  represents  the  0th 
iteration  of  the  solution  process.  The  terms  of  the  solution  may 
then  be  used  as  an  initial  guess  for  the  first  iteration  of  (23). 
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This  paper  is  the  second  in  a  two-part  study  of  identifying  mistuning  in  biaded  disks.  It 
presents  experimental  validation  of  a  new  method  of  mistuning  identification  based  on 
measurements  of  the  vibratory  response  of  the  ^stem  as  a  whole.  As  a  system-based 
method,  this  approach  is  particularly  suited  to  integrally  biaded  rotors,  whose  blades 
cannot  he  removed  for  individual  measurements.  The  method  is  based  on  a  recently 
developed  reduced-order  model  of  mistuning  called  the  fundamental  mistuning  model 
(FMM)  and  is  applicable  to  isolated  families  of  modes.  Two  versions  of  FMM  system 
identification  are  applied  to  the  experimental  data:  a  basic  version  . that  requires  some 
prior  knowledge  of  the  system*s  properties,  and  a  somewhat  more  complex  version  that 
determines  the  mistuning  completely  from  experimental  data.  [DOI:  10.1115/1. 16439 U] 


1  Introduction 

This  is  the  second  of  two  papers  on  identifying  mistiming  ip 
biaded  disks.  The  first  paper  [1]  reviews  the  literature,  derives  a 
new  theory  of  identifying  mistuning,  and  illustrates  its  applicabil¬ 
ity  using  numerical  examples.  Since  the  method  is  based  on  mea¬ 
surements  of  the  system  as  a  whole,  it  is  particularly  suited  to 
integrally  biaded  rotors  (IBRs).  Two  versions  of  system  identifi¬ 
cation  were  derived  in  the  first  paper:  a  basic  version  that  requires 
some  prior  knowledge  of  the  system’s  properties,  and  a  somewhat 
more  complex  version  that  determines  the  mistuning  completely 
from  experimental  data.  The  second  method  not  only  determines 
the  mistuning  in  the  IBR  but  also  determines  the  natural  frequen¬ 
cies  that  it  would  have  had  if  all  the  blade  and  disk  sectors  were 
identical,  i.e.,  the  “tuned”  system  frequencies.  In  this  second  pa¬ 
per,  we  apply  the  new  methods  to  a  modem  compressor  stage, 
identify  themistuning  and  the  tuned  system  frequencies,  and  show 
that  the  results  correlate  well  with  the  experimental  data  and  also 
with  independent  measurements  and  calculations  noade  by  engi¬ 
neers  at  Pratt  &  Whitney  Aircraft. 

The  literature  in  mistuning  is  extensive,  [2],  and  the  current 
research  is  put  in  context  in  Part  I,  [1].  However,  it  should  be 
emphasized  that  Judge  et  al.  introduced  the  concept  of  using  the 
system  modes  and  frequencies  to  determine  the  mistuning  in 
IBRs,  [3],  and  also  did  extensive  comparisons  with  experimental 
data.  The  key  differences  between  their  approach  and  that  reported 
here  are:  the  relative  simplicity  of  FMM  ID  and  the  consequence 
that  the  mistuning  in  the  system  can  be  determined  completely 
from  experimental  data;  that  FMM  ID  uses  a  blade/disk  sector 
measure  of  mistuning;  and  the  fact  that  their  method  is  not  inher¬ 
ently  limited  to  an  isolated  family  of  modes  as  is  FMM  ID. 

TTie  system  identification  methods  presented  here  will  be  ap¬ 
plied  to  two  IBRs  of  the  same  design.  The  IBRs  used  in  this  study 
were  designed  by  Prof.  S.  Fleeter  at  Purdue  University  in  coop¬ 
eration  with  Pratt  &  Whitney  Aircraft  to  reflect  the  aerodynamic 
and  structural  properties  of  a  modem  compressor.  The  work  re¬ 
ported  in  this  paper  is  part  of  a  major  research  initiative  on  mis¬ 
tuning,  friction  damping,  and  forced  response  that  is  sponsored  by 
the  U.S.  Air  Force,  the  U.S.  Navy,  NASA,  and  the  industrial  mem- 
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bers  of  the  GUIde  Consortium.  Consequently,  the  vibratory  re¬ 
sponse  of  the  compressor  stages  discussed  in  this  paper  will  be  the 
focus  of  mistuning  and  forced  response  studies  for  some  time  to 
come. 

This  paper  is  organized  as  follows.  Section  2  describes  the  two 
test  rotors  and  the  testing  procedure.  Then,  in  Section  3,  we 
present  the  FMM  ID  results.  The  determined  values  of  mistuning 
are  used  as  input  with  the  standard  FMM  prediction  code,  [4],  to 
predict  the  rotor’s  vibratory  response  to  a  traveling  wave  excita¬ 
tion  and  the  results  compared  with  experimental  data  in  Section  4. 
An  interesting  result  from  Section  3  is  that  the  mistuning  patterns 
in  both  IBRs  are  very  similar.  The  implications  of  this  result  are 
discussed  in  Section  5.  Finally,  the  key  attributes  of  the  method 
are  summarized  in  Conclusions. 

2  Test  Rotors  and  Procedures 

2.1  Test  Rotors.  To  investigate  how  well  FMM  ID  works 
with  real  experimental  data  from  actual  hardware,  we  applied  the 
methods  in  Part  1,  [1],  to  a  pair  of  transonic  compressors.  Fig.  1. 
The  two  rotors  are  designated  as  SN-1  and  SN-3. 

Our  industrial  partner  on  this  project,  Pratt  &  Whitney,  pro¬ 
vided  a  single  blade/disk  sector  fi^te  element  model  of  the  tuned 
compressor.  By  solving  this  model  with  free  boundary  conditions 
at  the  hub  and  various  cyclic  symmetric  boundary  conditions  on 
the  radial  boundaries  of  the  disk,  we  generated  a  nodal  diameter 
map  of  the  tuned  rotor.  Fig.  2.  The  free  boundary  conditions  at  the 
hub  represented  the  boundary  conditions  in  our  experiment:  an 
IBR  supported  by  a  soft  foam  pad  and  is  otherwise  unconstrained. 
Note  in  the  figure  that  each  of  the  first  two  families  of  modes  have 
isolated  frequencies.  These  correspond  to  first  bending  and  first 
torsion  modes,  respectively.  Since  FMM  ID  is  applicable  for  iso¬ 
lated  families  of  modes,  both  the  first  bending  and  first  torsion 
modes  are  suitable  candidates  for  our  identification  method 

2.2  Experimental  Procedures.  FMM  ID  requires  measure¬ 
ments  of  the  mistuned  rotor’s  system  modes  and  natural  frequen¬ 
cies.  By  system  mode,  we  mean  the  tip  displacement  of  each 
blade  as  a  function  of  angular  position.  These  inodes  were  ob¬ 
tained  using  a  standard  modal  analysis  approach:  measure  the 
biaded  disk’s  transfer  functions,  and  then  curve-fit  the  transfer 
functions  to  obtain  the  modes  and  natural  frequencies. 

Our  industrial  partners  performed  standard  transfer  function 
measurements.  The  rotor  was  placed  on  a  foam  pad  to  approxi- 
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Fig,  1  One  of  two  nominally  identical  test  compressors 


mate  a  free  boundary  condition.  Then,  they  excited  one  of  the 
blades  over  the  frequency  range  of  interest,  measured  the  response 
of  each  blade  with  a  laser  vibrometer,  and  determined  the  transfer 
functions  using  a  spectrum  analyzer.  A  typical  transfer  function  is 
shown  in  Fig.  3.  Note  that  due  to  the  high  modal  density,  it  was 
necessary  to  measure  the  response  with  a  very  high  frequency 
resolution.  This  process  was  repeated  for  both  compressors  over 
two  frequency  bands  in  order  to  capture  the  response  of  both  the 
first  bending  and  first  torsion  modes. 

We  then  used  the  commercially  available  MODENT  modal 
analysis  package  to  curve-fit  the  transfer  functions.  This  resulted 
in  measurements  of  the  mistuned  first  bending  and  torsion  modes 
of  each  rotor,  along  with  their  natur^  frequencies.  Because  the 
blade  that  was  excited  was  at  a  low  response  point  in  some  modes, 


Fig.  2  Natural  frequencies  of  compressor  with  no  mistuning 


Fig.  3  Representative  transfer  function  from  compressor  SN-1 


we  were  not  able  to  measure  two  or  three  of  the  modes  in  each 
family.  In  the  following  section,  we  use  these  measured  mistuned 
modes  and  natural  frequencies  to  demonstrate  the  applicability  of 
FMM  ID  to  actual  hardware. 


3  FMM  ED  Results 

The  measured  modes  and  frequencies  were  used  to  test  both 
forms  of  the  FMM  ID  method:  Basic  FMM  ID  and  the  completely 
experimental  Advanced  FMM  ID.  The  method  was  applied  to 
each  rotor,  for  both  the  first  bending  and  torsion  families  of 
modes.  The  tuned  frequencies  required  by  basic  FMM  ID  were  the 
same  as  those  depicted  in  Fig.  2. 

In  order  to  assess  the  accuracy  of  FMM  ID,  the  results  were 
compared  to  benchmark  data.  In  Section  3.1,  we  discuss  a  method 
for  obtaining  a  benchmark  measure  of  mistuning.  Then,  in  Sec¬ 
tions  3.2  and  3.3  we  present  the  results  of  FMM  ID  for  the  first 
bending  and  first  torsion  modes,  respectively. 

3.1  Benchmark  Measure  of  Mistuning.  In  order  to  assess 
the  accuracy  of  the  FMM  ID  method,  the  results  must  be  com¬ 
pared  to  benchmark  data.  However,  since  the  test  rotors  are  inte¬ 
grally  bladed,  their  mistuning  could  not  be  measured  directly. 
Therefore,  an  indirect  approach  was  used  to  obtain  the  benchmark 
mistuning.  Our  industrial  partners  carefully  measured  the  geom¬ 
etry  of  each  blade  on  the  two  rotors.  From  the  geometries,  they 
constmeted  finite  element  models  for  each  blade  and  calculated 
the  frequencies  that  it  would  have  if  it  were  clamped  at  its  root. 
Since  each  blade  had  a  slightly  different  geometry  it  also  had 
slightly  different  frequencies.  Thus,  the  variations  in  the  blade 
frequencies  caused  by  geometric  variations  were  determined.  This 
data  was  provided  to  CMU  and  we  put  it  in  a  form  that  could  be 
compared  with  the  values  identified  by  JMM  ID,  First,  we  calcu¬ 
lated  the  frequency  variations  as  a  fraction  of  the  mean  so  that  we 
knew  the  deviation  in  the  blade  frequencies.  These  in  turn  had  to 
be  related  to  the  sector  frequency  deviations  determined  by  FMM 
ID.  For  modes  with  most  of  their  strain  energy  in  the  blade,  sector 
frequency  deviations  can  be  obtained  from  blade  frequency  mis¬ 
tuning  by  simple  scaling,  i.e., 

Aw^=a(A<dj)  ,  (1) 

where  a  is  the  fraction  of  strain  energy  in  the  blade  for  the  average 
nodal  diameter  mode. 

Section  3.2  presents  the  results  for  the  family  of  bending 
modes,  and  Section  3.3  presents  the  results  from  the  torsion 
modes. 

3.2  FMM  ID  Results  for  Bending  Modes 

3.2.1  SN-1  Results.  The  measured  mistuned  modes  and 
natural  frequencies  for  the  compressor  SN-1  were  used  as  input  to 
both  versions  of  FMM  ID,  as  described  in  [1].  In  the  case  of  basic 
FMM  ID,  the  tuned  system  frequencies  of  the  first  bending  family 
from  Fig.  2  were  also  used  as  input. 

Figure  4  shows  the  sector  frequency  deviations  identified  by 
each  FMM  ID  method  along  with  the  benchmark  results.  In  both 
cases,  the  agreement  is  good.  This  implies  that  the  mistuning 
is  predominantly  caused  by  geometric  variations  and  that  the 
variations  are,  in  fact,  accurately  captured  by  Pratt  &  Whitney's 
process. 

In  order  to  make  these  comparisons  easier,  we  plotted  all  mis¬ 
tuning  in  Fig.  4  as  the  variation  from  a  zero  mean.  However,  it 
should  be  noted  that  this  rotor  had  a  mean  frequency  1.3%  higher 
than  that  of  the  tuned  finite  element  model.  This  DC  shift  was 
detected  by  Basic  FMM  ID  as  a  constant  amount  of  mistuning 
added  to  each  blade's  frequency.  However,  since  the  Advanced . 
FMM  ID  formulation  does  not  incorporate  the  tuned  finite  ele 
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Fig.  4  Comparison  of  mistuning  from  FMM  ID  with  benchmark 
results  for  SN-1 


ment  frequencies,  it  has  no  way  to  distinguish  between  a  mean 
shift  in  the  mistuning  and  a  corresponding  shift  in  the  tuned  sys¬ 
tem  frequencies.  Therefore,  in  Advanced  FMM  ID  we  define  mis¬ 
tuning  to  have  a  zero  mean,  arid  then  infer  a  corresponding  set  of 
tuned  frequencies. 

The  tuned  frequencies  identified  by  Advanced  FMM  ID  are 
compared  with  the  finite  element  values  in  Fig,  5.  Notice  that  the 
FMM  ID  frequencies  are  approximately  17  Hz  higher  than  the 
finite  element  values.  This  corresponds  to  a  1.3%  shift  in  the  mean 
of  the  tuned  system  frequencies  ^at  compensates  for  fact  that  the 
blade  mistuning  now  has  a  zero  mean.  To  facilitate  the  compari¬ 
son  of  the  finite  element  and  EMM  ID  results,  we  have  subtracted 
off  the  mean  shift,  and  then  plotted  the  results  as  circles  on  Fig.  5. 
Once  this  adjustment  is  made,  it  can  be  seen  that  the  distribution 
of  the  tuned  frequencies  determined  by  FMM  ID  agree  quite  well 
with  the  values  calculated  from  the  finite  element  model.  Clearly, 
the  finite  element  model  captures  the  same  variations  in  the  tuned 
system  frequencies  as  identified  by  Advanced  FMM  ID.  However, 
advanced  FMM  ID  identifies  the  fact  that  SN-1  had  slightly  higher 
average  frequencies  than  the  FEM  model — a  fact  that  could  be 
important  in  establishing  frequency  margins  for  the  stage. 


Fig.  5  Comparison  of  tuned  system  frequencies  from  FMM  ID 
and  FEM  for  SN-1 


Fig.  6  Comparison  of  mistuning  from  FMM  ID  with  benchmark 
results  for  SN-1 


Consider  the  sector  frequency  deviations  of  SN-1  shown  in  Fig. 
4.  Notice  that  the  mistuning  varies  from  blade-to-blade  in  a  regu¬ 
lar  pattern.  To  highlight  this  pattern,  we  will  reassign  the  blade 
numbering  so  that  blade  position  1  corresponds  to  the  high  fre¬ 
quency  blade,  Fig.  6.  Notice  that  when  plotted  in  the  new  num¬ 
bering  scheme,  the  mistuning  pattern  has  a  predominandy  de¬ 
creasing  trend,  with  a  jump  at  position  9.  This  trend  suggests  that 
the  mistuning  might  have  been  caused  by  tool  wear  during  the 
machining  process  and  that  an  adjustment  in  the  process  was 
made  after  half  of  the  18  blades  were  manufactured-  This  hypoth¬ 
esis  will  be  reexamined  after  reviewing  the  results  for  SN-3. 

3.2.2  SN-3  Results.  The  Basic  and  Advanced  FMM  ID 
methods  were  then  applied  in  a  similar  maimer  to  rotor  SN-3’s 
family  of  first  bending  modes.  The  identified  rinstuning  and  tuned 
system  frequencies  are  shown  in  Figs.  1  and  8.  For  comparison 
purposes,  we  have  again  plotted  the  mistuning  with  a  zero  mean, 
and  subtracted  a  corresponding  mean  shift  from  the  predicted 
tuned  system  frequencies.  The  agreement  is  also  good  for  rotor 
SN-3. 

In  Fig.  7,  we  numbered  the  blades  so  that  blade  1  corresponds 
to  the  high  frequency  sector.  Since  we  used  a  similar  numbering 
scheme  in  SN-1,  we  can  more  easily  compare  the  mistuning  in 


Fig.  7  Comparison  of  mistuning  from  FMM  ID  with  benchmark 
results  for  SN-3 
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Fig.  8  Comparison  of  tuned  system  frequencies  from  FMM  ID 
and  FEM  for  SN-3 


both  rotors.  It  is  interesting  to  note  that  the  mistuning  pattern  in 
Fig.  7  for  SN-3  is  quite  similar  to  that  of  Fig.  6  for  SN-1.  This 
result  proves  that  in  IBRs  mistuning  is  not  always  a  random  phe¬ 
nomenon.  The  implication  of  this  fact  on  the  predictability  of  the 
vibratory  response  is  discussed  in  Section  5. 


Fig.  9  Comparison  of  mistuning  from  FMM  ID  with  benchmark 
results  for  torsion  modes,  (a)  SN-1,  (d)  SN-3 
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Fig.  1 0  Comparison  of  tuned  system  frequencies  from  FMM  ID 
and  FEM  for  torsion  modes 


3.3  FMM  ID  Results  for  Torsion  Modes.  In  this  section 
we  will  examine  FMM  ID’s  ability  to  identify  mistuning  in  the 
first  torsion  modes.  Only  the  results  for  Advanced  FMM  ID  are 
presented  in  order  to  conserve  space. 

Advanced  FMM  ID  was  applied  to  each  test  rotor’s  family  of 
torsion  modes.  Figure  9  compares  the  mistuning  identified  by 
FMM  ID  with  the  values  inferred  by  our  industrial  partners  from 
their  geometric  measurements.  The  agreement  between  the  two 
methods  for  SN-1  is  quite  good,  while  the  agreement  for  SN-3  is 
remarkable.  In  Fig.  .9,  the  blades  are  numbered  in  the  same  order 
as  in  plots  6  and  7.  Thus,  the  mistuning  patterns  in  the  torsion 
modes  are  very  similar  to  those  observed  for  the  bending  modes, 
e.g.,  the  blades  with  the  highest  and  lowest  frequencies  are  the 
same  for  both  sets  of  modes.  This  suggests  that  the  mistuning  in 
these  systems  may  well  be  caused  by  relatively  uniform  thickness 
variations  in  the  blades  since  this  would  affect  the  frequencies  of 
both  types  of  modes  in  a  very  similar  manner. 

In  addition  to  identifying  the  mistuning  in  these  rotors,  Ad¬ 
vanced  FMM  ID  simultaneously  inferred  the  tuned  system  fre¬ 
quencies  of  the  system’s  torsion  modes,  as  shown  in  Fig.  10. 
Again,  the  agreement  is  good.  Thus,  FMM  ID  worked  well  on 
both  the  torsion  and  bending  modes  of  the  test  compressors. 


4  Forced  Response  Prediction 

In  the  previous  section,  we  found  that  both  forms  of  FMM  ID 
inferred  mistuning  that  agreed  very  well  with  benchmark  data.  In 
this  section,  the  identified  mistuning  will  be  used  to  predict  the 
forced  response  of  the  compressors  to  a  traveling  wave  excitation. 
The  results  will  be  compared  with  measurements  done  by  our 
industrial  partners. 

Pratt  &  Whitney  has  developed  an  experimental  capability  for 
simulating  traveling  wave  excitation  in  stationary  rotors.  Their 
technique  was  applied  to  SN-1  in  order  to  measure  its  first  bend¬ 
ing  family’s  response  to  a  3E  excitation.  We  then  predicted  the 
response  of  SN-1  with  the  methods  developed  here  and  in  [1,4]. 

The  issue  is:  do  the  28  parameters  (18  mistuned  frequencies  and 
10  tuned  system  frequencies)  identified  by  FMM  ID  from  one  set 
of  transfer  functions  determine  the  system  properties  sufficiently 
well  fiiat  we  can  accurately  predict  the  traveling  wave  response? 
To  make  the  prediction,  we  use  the  mistuning  and  tuned  system 
frequencies  from  Advanced  FMM  ID  as  input  to  the  FMM 
reduced-order  model.  FMM  calculates  the  system’s  mistuned 
modes  and  natural  frequencies.  Then,  we  use  modal  sununation  to 
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Fig.  11  Comparison  of  FMM  based  forced  response  with  ex¬ 
perimental  data,  (a)  FMM,  (b)  experiment 


calculate  the  response  to  a  3E  excitation.  The  modal  damping  used 
in  the  summation  was  calculated  from  the  half-power  bandwidth 
of  the  transfer  function  peaks. 

Figure  11  shows  the  comparison  of  the  benchmark  forced  re¬ 
sponse  results  with  that  predicted  by  FMM.  For  clarity,  only  the 
envelope  of  the  blade  response  is  shown.  Also,  the  plots  have  been 
normalized  so  that  the  maximum  response  is  equal  to  one.  In 
general,  the  two  curves  agree  reasonably  well.  In  order  to  obsen^e 
how  well  the  response  of  individual  blades  was  predicted  we  have 
also  compared  the  relative  responses  of  the  blades  at  two  resonant 
peaks,  the  peaks  labeled  ©  and  (2)  in  Fig,  11.  The  relative  ampli¬ 
tude  of  each  blade  as  determined  by  both  methods  is  plotted  for 
both  resonant  peaks  in  Fig.  12.  The  agreement  is  also  reasonably 
good.  Thus,  the  FMM  based  method  not  only  captured  the  overall 
shape  of  the  response,  but  also  determined  the  relative  amplitudes 
of  the  blades  at  the  various  resonances. 

5  Cause  and  Implications  of  Repeated  Mistuning  Pat¬ 
tern 

In  the  literature,  the  mistuning  in  bladed  disks  is  generally  con¬ 
sidered  to  be  a  random  phenomenon.  However,  in  Section  3  we 
saw  that  both  test  rotors  have  very  similar  mistuning  patterns  that 
are  far  from  random.  If  such  repeated  mistuning  matters  are  found 
to  be  common  among  IBRs,  it  will  have  broad  implications  on  the 
predictability  of  these  systems.  In  Section  5.1  we  discuss  the 


Blade  Position 


Fig.  12  Relative  blade  amplitudes  at  forced  response  reso¬ 
nance,  (a)  Resonance  1,  (b)  Resonance  2 


cause  of  the  repeated  mistuning.  Then  in  Section  5.2,  we  examine 
its  implications  on  forced  response  predictability  for  a  larger 
sample. 

5.1  Possible  Cause  of  Repeated  Mistuning.  The  similarity 
between  the  mistuning  patterns  identified  in  SN-1  and  SN-3  is 
highly  suggestive  that  the  mistuning  was  caused  by  a  consistent 
manufacturing  effect.  In  addition,  we  observed  that  the  mistuning 
in  the  torsion  modes  follows  the  same  trends  as  in  the  bending 
modes.  Thus,  the  dominant  form  of  mistuning  is  most  likely 
caused  by  relatively  uniform  blade-to-blade  thickness  variations. 
One  plausible  explanation  for  the  observed  patterns  is  tool  wear. 
Suppose  that  the  blades  were  machined  in  descending  order  from 
blade  18  to  blade  1.  Then,  due  to  tool  wear,  each  subsequent  blade 
will  be  slightly  larger  than  the  previous  one.  This  effect  would 
cause  the  sector  frequencies  to  inonotonically  increase  around  the 
wheel.  With  the  exception  of  the  frequency  jump  observed  at 
blade  9,  this  behavior  matches  the  observed  mistuning.  Since 
there  are  18  blades,  the  discontinuity  at  blade  9  could  well  be  the 
result  of  a  tool  adjustment  made  hdfway  through  the  machining 
process. 

5.2  Implications  of  Repeated  Mistuning.  The  repeating 
mistuning  patterns  caused  by  such  machining  effects  can  signifi¬ 
cantly  increase  our  abOity  to  accurately  predict  the  response  of  the 
fleet  through  probabilistic  methods.  For  example^  consider  an  en- 
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Fig.  13  CDF  of  peak  blade  amplitude  for  a  nominally  tuned 
and  nominally  mistuned  compressor 


tire  fleet  of  the  transonic  compressors  studied  in  this  paper.  If  we 
were  to  inconectly  assume  that  the  mistuning  in  these  rotors  was 
completely  random,  then  we  would  estimate  that  the  sector  fre¬ 
quency  deviation  of  each  sector  has  a  mean  of  zero  and  a  standard 
deviation  of  about  2%.  Assuming  these  variations,  we  used  FMM 
to  perform  10,000  Monte  Carlo  simulations  to  represent  how  a 
fleet  of  engines  would  respond  to  a  3E  excitation.  We  used  the 
data  to  compute,  the  cumulative  distribution  function  (CDF)  of  the 
maximum  blade  amplitude  on  each  compressor.  The  CDF  of  a 
fleet  of  engines  with  random  mistuning  have  a  standard  deviation 
of  2%  is  shown  as  the  dashed  line  in  Fig.  13.  Notice  that  the 
maximum  amplitude  varies  widely  across  the  fleet,  ranging  in 
magnification  from  1.1  to  2.5. 

However,  these  rotors  are  in  fact  nominally  mistuned  with  a 
small  random  variation  about  the  nominal  pattern.  Since  the  ran¬ 
dom  variation  is  much  smaller  than  that  considered  above,  the 
fleet’s  response  is  actually  far  more  predictable.  To  illustrate  this 
point,  we  approximated  the  nominal  mistuning  pattern  as  the 
mean  of  the  patterns  measured  for  the  two  test  rotors.  Based  on 
this  pattern,  we  found  that  the  sector  frequency  deviations  differed 
from  the  nominal  values  with  a  standard  deviation  of  only  0.2%, 
as  shown  in  Fig.  14,  Making  use  of  the  fact  that  the  rotors  are 
nominally  mistuned,  we  repeated  the  Monte  Carlo  simulations. 
We  then  computed  the.  CDF  of  the  maximum  amplitude  on  each 
rotor.  The  results  are  plotted  as  the  solid  line  on  Fig.  13,  Notice 
that  by  accounting  for  nominal  mistuning,  the  range  of  maximum 
amplitudes  is  significantly  reduced.  Thus,  if  we  can  measure  and 
make  use  of  nominal  mistuning  when  it  occurs  then  the  fleet’s 
behavior  will  be  far  more  predictable, 

6  Conclusions 

A  new  method  of  identifying  mistuning  in  bladed  disks  is  dem¬ 
onstrated  using  actual  hardware.  The  method  is  called  FMM  ID 
because  it  is  based  on  the  fundamental  mistuning  model  (FMM), 
[4].  To  test  the  FMM  ID  approach,  we  used  it  to  identify  the 
mistuning  in  a  pair  of  modem  transonic  compressors  and  com¬ 
pared  the  results  with  values  that  were  determined  by  Pratt  & 
Whitney  using  a  completely  independent  method. 

FMM  ID  uses  measurements  of  the  system  mode  shapes  and 
natural  frequencies  to  infer  the  rotor’s  mistuning.  The  key  concept 
behind  FNM  ID  is  that  the  high  sensitivity  of  system  modes  to 
small  variations  in  mistuning  causes  measurements  of  those 
modes  themselves  to  be  an  accurate  basis  for  mistuning  identifi¬ 
cation.  Since  FMM  ID  does  not  require  individual  blade  measure- 


Fig.  14  Mean  and  standard  deviation  of  each  sector’s  mistun- 
ing  for  a  nominally  mistuned  compressor 


ments,  it  is  particularly  suited  to  integrally  bladed  rotors.  The 
method  is  only  applicable  to  isolated  families  of  modes.  We  have 
developed  two  forms  of  FMM  ID:  Basic  FMM  ID,  and  an  ad¬ 
vanced  version  that  also  identifies  the  frequencies  of  the  tuned 
system.  - 

We  applied  both  forms  of  FMM  ID  to  the  two  test  rotors.  To 
provide  benchmark  values  of  mistuning,  our  industrial  partners  in 
this  research,  Pratt  &  Whitney,  also  identified  the  mistuning  in  the 
rotors  by  measuring  the  geometry  of  each  blade  and  then  deter¬ 
mining  its  natural  frequencies  through  a  finite  element  analysis. 
The  results  from  both  methods  of  identifying  mistuning  agreed 
veiy  well.  In  addition  to  identifying  mistuning,  the  advanced  form 
of  FMM  ID  can  also  infer  the  tuned  system  frequencies  of  the 
rotor.  Except  for  a  slight  shift  in  their  mean  value,  the  identified 
tuned  frequencies  agreed  weU  with  finite  element  values.  In  prac¬ 
tice,  the  ability  of  FMM  ID  to  identify  the  tuned  frequencies  of 
the  system  could  provide  a  useful  method  of  confirming  that  the 
manufacturing  process  resulted  in  an  IBR  that  is  consistent  with 
the  original  design. 

We  observed  that  the  mistuning  patterns  in  the  two  test  rotors 
were  very  similar.  This  suggests  that  the  mistuning  was  caused  by 
a  repeating  feature  of  the  manufacturing  process,  perhaps  tool 
wear.  The  cause  of  the  repeating  mistuning  pattern  will  be  the 
subject  of  further  investigation.  If  the  mistuning  in  the  two  IBRs 
that  we  have  tested  are,  in  fact,  representative  of  the  manufactur¬ 
ing  process,  then  a  larger  sample  of  rotors  would  exhibit  similar 
behavior,  i.e,,  they  would  have  significant  levels  of  nominal  mis¬ 
tuning  with  small  levels  of  random  mistuning  superimposed.  In 
this  case,  the  forced  response  of  a  fleet  of  these  compressors 
would  be  far  more  predictable  than  might  have  been  previously 
foreseen.  If  the  manufacturing  process  can  be  understood  and  con¬ 
trolled,  then  we  may  be  able  to  use  optimization  techniques  such 
as  those  proposed  in  [5,6]  to  manufacture  IBRs  that  have  low, 
robust  response. 

The  FMM  ID  method  provides  a  simple  approach  for  accu¬ 
rately  identifying  mistuning  in  integrally  bladed  rotors  for  isolated 
families  of  modes.  la  fact,  because  of  its  simplicity  one  version  of 
the  method  requires  only  experimental  data  to  determine  the  key 
parameters  that  characterize  its  response.  Once  determined,  the 
mistuning  can  be  used  with  the  FMM  reduced-order  model,  [4],  to 
predict  how  much  the  mistuning  will  increase  the  stage’s  forced 
response  to  a  traveling  wave,  engine  order  excitation.  This  tech¬ 
nology  is  useful  since  it  will  allow  test  engineers  to  determine 
how  vibratory  response  of  a  specific  IBR  that  is  tested  in  a 
spin  pit,  rig  test  or  engine  relates  to  the  vibratory  response  of  the 
population  as  a  whole.  Ultimately,  this  technology  should  allow  us 
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to  reduce  the  number  of  tests  required  to  characterize  the  vibratory 
response  of  components  and  still  have  more  durable  engines. 
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ABSTRACT 

A  new  approach  to  modal  analysis  is  presented  that  allows 
the  modes  and  natural  frequencies  of  a  mistimed  bladed  disk  to 
be  determined  from  its  response  to  a  traveling  wave  excitation. 
The  resulting  modes  and  natural  frequencies  are  then  used  as 
input  to  a  system  identification  method  to  determine  the  bladed 
disk's  mistuning  while  it  is  rotating.  This  capability  is  useful 
since  it  provides  a  basis  for  determining  blade  frequencies 
under  engine  operational  conditions  and  could  help  monitor  the 
health  of  the  engine. 

INTRODUCTION 

Bladed  disks  used  in  turbine  engines  are  nominally 
designed  to  be  cyclically  symmetric.  If  this  were  the  case,  then 
all  blades  would  respond  with  the  same  amplitude  when  excited 
by  a  traveling  wave.  However,  in  practice,  the  resonant 
amplitudes  of  the  blades  are  very  sensitive  to  small  changes  in 
their  properties.  These  variations  are  referred  to  as  mistuning, 
and  may  result  from  the  manufacturing  process  or  wear. 
Mistuning  causes  some  blades  to  have  a  significantly  higher 
vibratory  response  and  may  cause  them  to  fail  from  high  cycle 
fatigue.  Thus,  mistuning  is  a  primary  source  of  uncertainty  in  a 
bladed  disk  system's  response.  Srinivasan  provides  a  review  of 
this  topic  in  [1]. 

Two  types  of  tools  have  recently  been  developed  to  help 
predict  the  response  of  mis  tuned  systems  [2-5]:  reduced  order 
models  (ROMs)  and  system  identification  methods.  The  ROMs 
are  very  efficient  methods  for  accurately  predicting  the  forced 
response  of  a  bladed  disk.  They  provide  structural  fidelity 
comparable  to  a  finite  element  analysis  of  a  full  mistuned 


bladed  disk  with  the  numerical  efficiency  similar  to  that  of  a 
simple  mass  spring  model.  Thus,  they  can  be  used  as  the  basis 
of  Monte  Carlo  probabilistic  analyses  and  optimization 
programs.  In  addition,  they  reduce  the  number  of  parameters 
that  characterize  mistuning  to  a  manageable  level  and  become 
the  basis  of  system  identification  methods. 

System  identification  methods  can  determine  the  mistuning 
of  a  specific  bladed  disk.  Two  methods  have  been  reported  in 
the  literature  for  determining  mistuning  from  the  vibratory 
response  of  the  coupled  bladed  disk  system.  The  methods  are 
based  on  different  ROMs  of  the  mistuned  system.  The  first  is 
derived  from  the  reduced  order  model  REDUCE  [3]  and  the 
second  from  FMM  [5].  For  the  purpose  of  system 
identification,  the  equations  in  the  ROMs  are  reformulated  so 
that  they  can  be  used  to  solve  the  inverse  problem,  i.e.  given  the 
mistuned  modes  and  natural  frequencies  of  the  system 
determine  the  amount  of  mistiming  in  each  blade/disk  sector. 

The  first  method  to  use  system  modes  to  determine 
mistuning  was  based  on  REDUCE.  Since  REDUCE  was 
derived  using  component  mode  synthesis,  it  requires  a 
significant  amount  of  preliminary  finite  element  analysis  before 
it  can  be  used  for  system  identification  [6].  The  second  method 
was  based  on  the  Fimdamental  Mistuning  Model  (FMM)  and  is 
called  FMM  ID  [7,  8].  The  advantage  of  this  method  is  that  it  is 
much  simpler  and,  as  a  result,  the  mistuning  in  the  system  can 
be  identified  without  any  preliminary  modeling  of  the  system, 
i.e.  the  method  is  completely  experimental.  Its  disadvantage  is 
that  it  only  works  if  the  family  of  modes  under  consideration 
has  frequencies  that  are  relatively  isolated.  In  both  methods  of 
system  identification,  the  mistuned  modes  and  natural 
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frequencies  of  the  bladed  disk  are  used  as  input  data  to 
determine  the  mistuning  in  the  system. 

These  recent  system  identification  methods  overcome  a 
significant  limitation  of  the  traditional  approach  for  determining 
a  bladed  disk’s  mistuning.  Traditionally,  mistuning  in  rotors 
with  attachable  blades  is  measured  by  mounting  each  blade  in  a 
broach  block  and  measuring  its  natural  firequency.  The 
difference  of  each  blade’s  frequency  from  the  mean  value  is 
then  taken  as  a  measure  of  its  mistuning.  However,  this  method 
cannot  be  applied  to  integrally  bladed  rotors  (IBRs)  whose 
blades  cannot  be  removed  for  individual  testing.  In  contrast, 
the  recent  system  identification  techniques  rely  on 
measurements  of  the  bladed  disk  system  as  a  whole,  and  are 
thus  well  suited  to  IBRs. 

These  modem  methods  are  also  potentially  valuable  tools 
for  determining  the  mistuning  in  conventional  bladed  disks. 
Even  when  applied  to  bladed  disks  with  conventionally  attached 
blades,  the  traditional  broach  block  method  of  mistuning 
identification  is  limited.  In  particular,  it  does  not  take  into 
account  the  fact  that  the  mistuning  measured  in  the  broach 
block  may  be  significantly  different  from  the  mistuning  that 
occurs  when  the  blades  are  mounted  on  the  disk.  This  variation 
can  arise  because  each  blade’s  frequency  is  dependent  on  the 
contact  conditions  at  the  attachment.  In  the  engine,  the 
attachment  is  loaded  by  centrifugal  force  from  the  blade  which 
provides  a  different  contact  condition  than  the  clamping  action 
used  in  broach  block  tests.  This  difference  is  accentuated  in 
multi-tooth  attachments  since  different  teeth  may  come  in 
contact  depending  on  how  the  attachment  load  is  applied.  In 
addition,  the  contact  in  multi-tooth  attachments  may  be 
sensitive  to  manufacturing  variations  and,  consequently,  vary 
from  one  location  to  the  next  on  the  disk.  To  address  these 
issues,  a  method  of  system  identification  is  needed  that  can  be 
used  to  directly  determine  mistuning  while  the  stage  is  rotating. 
The  method  will  need  to  identify  mistuning  from  the  response 
of  the  entire  system  since  the  blades  are  inherently  coupled 
under  rotating  conditions.  This  paper  presents  such  a  method. 

The  method  presented  in  this  paper  provides  an  approach 
for  extracting  the  mistuned  modes  and  natural  frequencies  of  the 
bladed  disk  under  rotating  conditions  from  its  response  to 
naturally  occurring,  engine  order  excitations.  The  method  is 
not  a  new  modal  analysis  method  but  rather  a  coordinate 
transformation  that  makes  traveling  wave  response  data 
compatible  with  the  existing,  proven  modal  analysis  algorithms. 
Once  the  modes  and  natural  frequencies  are  known  they  can  be 
used  as  input  to  either  method  of  system  identification.  This 
paper  will  use  FMM  ID  to  demonstrate  the  process. 

This  paper  is  organized  as  follows.  Section  2  presents  the 
theoretical  basis  of  die  mode  extraction  approach  for  use  with 
traveling  wave  excitations.  Then  in  Section  3,  we  present  two 
experimental  test  cases  that  illustrate  the  method.  In  each  case, 
the  extracted  modes  are  used  with  FMM  ID  to  determine  the 


rotor’s  mistuning,  which  is  then  compared  with  benchmark 
values.  Finally,  the  results  are  summarized  in  Conclusions. 

2.  THEORY 

Both  of  the  mistuning  identification  methods  cited  require 
the  mistuned  modes  and  natural  frequencies  of  the  bladed  disk 
as  input.  Under  stationary  conditions,  they  can  be  determined 
by  measuring  the  transfer  functions  of  the  system  and  using 
standard  modal  analysis  procedures.  One  way  of  measuring  the 
transfer  functions  is  to  excite  a  single  point  with  a  known 
excitation  and  measure  the  frequency  response  of  all  of  the 
other  points  that  define  the  system.  However,  when  the  bladed 
disk  is  subjected  to  an  engine  order  excitation  all  of  the  blades 
are  simultaneously  excited  and  it  is  not  clear  how  the  resulting 
vibratoiy  response  can  be  related  to  the  transfer  functions 
typically  used  for  modal  identification.  It  is  shown  here  that  if 
the  blade  frequency  response  data  is  transformed  in  a  particular 
manner  then  the  traveling  wave  excitation  constitutes  a  point 
excitation  in  the  transform  space  and  that  standard  modal 
analysis  techniques  can  then  be  used  to  extract  the  transformed 
modes.  Once  the  transformed  modes  are  determined  the 
physical  modes  of  the  system  can  be  calculated  from  an  inverse 
transformation. 

Section  2.1  describes  the  traditional  Single-Input-Single- 
Output  (SISO)  modal  analysis  method  and  identifies  its 
limitations  with  respect  to  multi-point  excitation  data  such  as 
the  response  to  a  traveling  wave.  In  section  2.2  a  general 
coordinate  transformation  is  presented  that  can  be  used  to 
express  multi-point  excitation  response  data  in  a  form  that  is 
compatible  with  traditional  SISO  modal  analysis  techniques. 
Then  in  section  2.3  we  describe  how  the  coordinate 
transformation  may  be  simplified  for  the  case  of  traveling  wave 
response  data. 

2.1  Traditional  Modal  Analysis 

Standard  modal  analysis  techniques  are  based  on 
measurements  of  a  structure’s  frequency  response  functions 
(FRFs).  These  frequency  response  functions  are  then 
assembled  as  a  frequency  dependant  matrix,  H(<2>) ,  in  which  the 
element  H^  j{a))  corresponds  to  the  response  of  point  i  to  the 
excitation  of  point  j  [9].  Modal  analysis  methods  require  that 
one  row  or  column  of  this  frequency  response  matrix  be 
measured.  In  the  cases  considered  in  this  paper,  the  mistuned 
modes  correspond  to  a  single  isolated  family  of  modes.  For 
example,  the  lower  frequency  modes  such  as  first  bending  and 
first  torsion  families  often  have  frequencies  that  are  relatively 
isolated.  When  this  is  the  case  the  “modes”  of  interest  are 
defined  in  terms  of  how  the  blade  displacements  vary  from  one 
blade  to  the  next  around  the  wheel  and  can  be  characterized  by 
the  response  of  one  point  per  blade.  Thus,  the  standard  modal 
analysis  experiment  may  be  performed  in  one  of  two  ways  when 
measuring  the  mistuned  modes  of  a  bladed  disk.  First,  the 
structure’s  frequency  response  may  be  measured  at  one  point  on 
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each  blade,  while  :t  is  excited  at  only  one  blade.  This  would 
result  in  the  measurement  of  a  single  column  of  R(a)) . 
Alternatively,  a  row  of  R((o)  may  be  obtained  by  measuring 
the  structure’s  response  at  only  one  blade  and  exciting  the 
system  at  each  blade  in  turn.  Notice  that  in  either  of  these 
acceptable  test  configurations,  the  structure  is  excited  at  only 
one  point  at  a  time.  However,  in  a  traveling  wave  excitation,  all 
blades  are  excited  simultaneously.  Thus,  the  response  of 
systems  subjected  to  such  multi-point  excitations  cannot  be 
directly  analyzed  by  standard  SISO  modal  analysis  methods. 

2.2  General  Multi-Point  Excitation  Analysis 

Since  a  traveling  wave  excitation  is  not  directly  compatible 
with  standard  SISO  modal  analysis  methods,  a  different 
approach  that  allows  for  multi-point  excitations  is  needed.  A 
significant  amount  of  research  has  been  performed  on  the  topic 
of  multiple-input  data  analysis.  However,  the  work  has  tended 
to  focus  on  the  case  of  random  excitations,  in  which  each 
excitation  point  drives  the  structure  with  a  different  set  of 
randomly  selected  frequencies  [10-12],  In  contrast,  a  traveling 
wave  excites  each  measurement  point  with  the  same  frequency 
at  any.  given  time.  The  method  presented  here  is  applicable  to 
any  multi-input  system,  in  which  the  frequency  profile  is 
consistent  from  one  excitation  point  to  the  next;  however,  the 
amplitude  and  phase  of  the  excitation  sources  are  free  to  vary 
spatially.  Suitable  excitation  forms  include  traveling  waves, 
acoustic  pressure  fields,  and  even  shakers  when  appropriately 
driven. 

In  typical  applications,  the  ij  element  of  the  frequency 
response  matrix  corresponds  to  the  response  of  point  i  to 

the  excitation  of  point  y.  However,  in  order  to  analyze 
frequency  response  data  from  a  multi-point  excitation,  we  must 
view  R^jico)  in  a  more  general  fashion.  In  a  more  general 
sense,  the  iJ  element  describes  the  response  of  the  coordinate 
to  an  excitation  at  the  coordinate.  Although  these 
coordinates  are  typically  taken  to  be  the  displacement  at  an 
individual  measurement  point,  this  need  not  be  the  case. 

The  structure’s  excitation  and  response  can  instead  be 
transformed  into  a  different  coordinate  system.  For  example,  an 
N  degree-of-freedom  coordinate  system  can  be  defined  by  a  set 
of  N  orthogonal  basis  vectors  which  span  the  space.  In  this 
representation,  each  basis  vector  is  a  coordinate.  The  key  to 
performing  modal  analysis  on  multi-point  excitation  data  is  to 
select  a  coordinate  system  in  which  the  excitation  is  described 
by  just  one  basis  vector.  Thus,  within  this  newly  defined  modal 
analysis  coordinate  system,  the  structure  is  subjected  to  only  a 
single  coordinate  excitation.  Therefore,  when  our  response 
measurements  are  expressed  in  this  same  domain,  they  represent 
a  single  column  of  the  FRF  matrix,  and  can  be  analyzed  by 
standard  SISO  modal  analysis  techniques.  The  following 
section  describes  how  this  approach  may  be  applied  to  traveling 
wave  excitations. 


2.3  Traveling  Wave  Modal  Analysis 

Consider  an  A/^bladed  disk  subjected  to  a  traveling  wave 
excitation.  It  is  assumed  that  the  amplitude  and  phase  of  each 
blade’s  response  is  measured  as  a  function  of  excitation 
frequency.  In  practice,  these  measurements  could  be  made 
under  rotating  conditions  with  a  Non-intrusive  Stress 
Measurement  System  (NSMS),  whereas  a  laser  vibrometer 
could  be  used  in  a  stationary  bench  test.  For  simplicity,  we  will 
only  consider  one  measurement  point  per  blade. 

It  is  assumed  that  the  blades  are  excited  harmonicly  by  the 
force  where  the  vector  f  describes  the  spatial 

distribution  of  the  excitation  force.  Similarly,  the  response  of 
each  measurement  point  is  given  by  h(a))e^‘"  .  The  components 
of  f  and  h  are  complex  since  they  contain  phase  as  well  as 
magnitude  information.  It  is  this  excitation  and  response  data 
from  which  we  wish  to  extract  modes  shapes  and  natural 
frequencies.  However,  in  order  for  this  data  to  be  compatible 
with  standard  SISO  modal  analysis  methods,  it  must  first  be 
transformed  to  an  appropriate  modal  analysis  coordinate 
system. 


As  indicated  in  Section  2.2,  an  appropriate  coordinate 
system  that  would  allow  this  to  occur  is  one  in  which  the  spatial 
distribution  of  the  force,  f ,  is  itself  a  basis  vector.  The  spatial 
distribution  of  a  traveling  wave  excitation  has  the  form:^ 


e 
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where  E  is  the  engine  order  of  the  excitation.  Therefore,  a 
coordinate  system  whose  basis  vectors  are  the  N  possible  values 
of  f ,  corresponding  to  all  N  distinct  engine  order  excitations,  0 
through  A^-1  may  be  used  as  a  basis.  The  basis  vectors  are 
complete  and  orthogonal. 


The  vectors  f  and  h  are  transformed  into  this  modal 
analysis  coordinate  system  by  expressing  them  as  a  sum  of  the 
basis  vectors.  Denoting  our  basis  vectors  as  the  set 
{bo,  bj , . . . ,  b^Y-i }  j  litis  summation  takes  the  form, 


N-\ 

njstO 

(2a) 

1 

m-Q 

(2b) 

‘Note,  that  in  order  to  simplify  the  concepts,  we  have  included  only  the 
phase  difference  that  occurs  from  one  blade  to  the  next.  In  the  case  of  higher 
frequency  applications  it  would  be  necessary  to  also  include  the  spatial 
variation  of  the  force  over  the  airfoil  if  more  than  one  family  of  modes  interact. 
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where  the  coefficients  and  h{co)^  describe  the  value  of  the 
w*  coordinate  in  the  modal  analysis  domain.  To  identify  the 
values  of  these  coefficients,  we  could  use  orthogonality;  this  is  a 
general  approach  that  is  applicable  for  any  orthogonal 
coordinate  system.  However,  for  the  case  of  traveling  wave 
excitations,  the  coordinate  transformation  may  be  simplified. 

Consider  the  element  of  the  vectors  in  (2).  For 
convenience,  let  all  vector  indices  run  from  0  to  iV-1.  Thus, 
these  elements  may  be  expressed  as, 

=  (3a) 

f»T=0 

(3b) 

/7r=0 

where  the  exponential  term  is  the  n'*  component  of  the  basis 
vector  b^.  Observe  that  (3)  is  the  inverse  discrete  Fourier 
Transform  (DFT‘^)  of  / .  This  relation  allows  us  to  state  the 
transformation  between  physical  coordinates  and  the  modal 
analysis  domain  in  the  simpler  form, 

f  =  DFT'‘{f}  (4a) 

h=:DFT"'{h}  (4b) 

and  conversely, 

f  =  DFT{f}  (5a) 

h  =  DFT{h}  (5b) 

where  DFT  is  the  discrete  Fourier  Transform  of  the  vector. 

By  applying  equation  (5),  the  force  and  response  vectors 
are  transformed  to  the  modal  analysis  coordinate  system.  Due 
to  our  selection  of  basis  vectors,  the  resulting  vector  f  will 
contain  only  one  nonzero  term  that  corresponds  to  the  engine 
order  of  the  excitation,  i.e.  a  5E  excitation  will  produce  a 
nonzero  term  in  element  5  of  f  .  This  indicates  that  within  the 
modal  analysis  domain,  we  have  excited  only  the  coordinate. 
Therefore,  h(fl>)  represents  column  E  of  the  FRF  matrix. 

The  transformed  response  data,  h(^),  may  now  be 
analyzed  using  standard  SISO  modal  analysis  algorithms.  The 
resulting  modes  will  also  be  in  the  modal  analysis  coordinate 
system,  and  must  be  converted  back  to  physical  coordinates 
though  an  inverse  discrete  Fourier  Transform,  (4).  These 
identified  modes  and  natural  frequencies  may  in  turn  be  used  as 
inputs  to  FMM  ID  to  determine  the  mistuning  of  a  bladed  disk 
from  its  response  to  an  engine  order  excitation. 

There  are  two  further  details  of  this  method  that  should  be 
noted.  First,  for  the  purpose  of  notation  convenience  we  have 


numbered  the  indices  of  all  matrices  and  vectors  from  0  to  iV'-l. 
However,  most  modal  analysis  packages  use  a  numbering 
convention  that  starts  at  1.  Therefore,  an  E!^  coordinate 
excitation  in  our  notation,  corresponds  to  an  (£'+l)^*  coordinate 
excitation  in  the  standard  convention.  This  must  be  taken  into 
account  when  specifying  the  “excitation  point”  in  the  modal 
analysis  software.  Second,  the  coordinate  transformation  used 
in  this  method  is  based  on  a  set  of  complex  basis  vectors.  Since 
the  modes  are  extracted  in  the  modal  analysis  domain  they  will 
be  highly  complex,  even  for  lightly  dmnped  systems.  Thus  it  is 
necessary  to  use  a  modal  analysis  package  that  can  properly 
handle  highly  complex  mode  shapes.  Not  all  commercial 
algorithms  are  suitable.  In  our  experience,  we  found  that  the 
MODENT  Suite  by  ICATS  [13]  works  quite  well  for  this 
purpose. 

The  next  section  presents  two  experimental  examples  of  the 
method  and  shows  how  it  can  be  used  for  system  identification. 

3,  EXPERIMENTAL  TEST  CASES 

This  section  presents  two  experimental  test  cases  of  the 
traveling  wave  system  identification  technique.  In  the  first 
example,  we  excited  an  integrally  bladed  fan  (IBR)  with  a 
traveling  wave  while  it  was  in  a  stationary  configuration. 
Because  the  IBR  was  stationary,  we  were  able  to  make  very 
accurate  response  measurements  using  a  laser  vibrometer. 
Thus,  this  example  serves  as  a  benchmark  test  of  the  traveling 
wave  identification  theory.  Then,  in  the  second  example  we 
explore  the  method’s  effectiveness  on  a  rotor  that  is  excited  in  a 
spin  pit  under  rotating  conditions.  The  amplitude  and  phase  of 
the  response  are  measured  using  an  NSMS  system;  NSMS  is  a 
non-contacting  measurement  method  which  is  commonly  used 
in  the  gas  turbine  industry  for  rotating  tests.  The  purpose  of  this 
example  is  to  assess  if  NSMS  technology  is  sufficiently 
accurate  that  it  can  be  used  with  our  traveling  wave  system 
identification  technique  to  determine  the  IBR’s  mistuning  from 
its  engine  order  response. 

3.1  Stationary  Benchmark 

Consider  the  integrally  bladed  fan  shown  in  Fig.  1.  This 
fan  was  tested  using  the  traveling  wave  excitation  system  at 
Wright  Patterson  Air  Force  Base’s  Turbine  Engine  Fatigue 
Facility  [14].  Since  the  facility’s  test  system  uses  an  array  of 
phased  electromagnets  to  generate  a  traveling  wave  excitation, 
the  bladed  disk  remains  stationary  during  the  test.  This 
configuration  is  ideal  for  benchmark  mistuning  studies  since  it 
allows  us  to  use  laser  vibrometry  to  obtain  very  accurate 
measurements  of  the  rotor’s  response. 

This  example  serves  two  purposes.  First,  it  validates  the 
traveling  wave  modal  analysis  method  presented  in  this  paper, 
and  second  it  demonstrates  that  the  resulting  modes  and  natural 
frequencies  may  be  used  to  determine  the  bladed  disk’s 
mistuning.  The  experiment  was  performed  with  the  fan  placed 
on  a  rubber  mat  to  approximate  a  free  boundary  condition. 
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Figure  1 :  Fan  tested  for  stationary  benchmark 


Figure  2:  Masses  used  to  mistune  fan 


First,  we  intentionally  mistimed  the  IBR  by  fixing  a 
different  mass  to  the  leading  edge  tip  of  each  blade  with  wax, 
Fig.  2.  The  masses  ranged  between  0  and  7g,  and  were  selected 
randomly.  Then,  to  obtain  a  benchmark  measure  of  the 
mistuned  fan’s  mode  shapes,  we  performed  a  standard  SISO 
modal  analysis  test.  Specifically,  we  used  a  single 
electromagnet  to  excite  blade  1  over  the  frequency  range  of  the 
first  bending  modes  while  we  measured  the  response  at  all  16 
blades  with  a  Scanning  Laser  Doppler  X^brometer  (SLDV). 
The  modes  were  then  extracted  from  the  measured  FRFs  using 
the  commercially  available  MODENT  modal  analysis  package. 
Figure  3  shows  a  portion  of  a  representative  FRF  of  this  system. 
Notice  that  the  modes  are  very  densely  packed,  however  due  to 
extremely  light  damping  we  were  able  to  extract  all  of  the 
modes  within  the  desired  frequency  range. 

Next,  to  validate  the  traveling  wave  modal  analysis  method 
the  fan  was  excited  using  a  5*  engine  order  traveling  wave 
excitation.  Again,  the  response  of  each  blade  was  measured 


using  the  SLDV.  The  blade  responses  to  the  traveling  wave 
excitation  were  transformed  using  equation  (5)  and  then 
analyzed  with  MODENT  to  extract  the  transformed  modes. 
Since  MODENT  numbers  its  coordinates  starting  at  1  (OE),  a 
5E  excitation  corresponds  to  the  excitation  of  coordinate  6. 
Therefore,  in  the  mode  extraction  process,  we  specified  that  the 
excitation  was  applied  at  the  6*  coordinate.  Lastly,  equation  (4) 
was  used  to  transform  the  resulting  modes  back  to  physical 
coordinates. 


Excitatton  Frequency  (Hz) 

Figure  3:  Representative  FRF 

The  modes  measured  through  the  traveling  wave  test  were 
then  compared  with  those  from  our  benchmark  analysis.  If  this 
system  were  tuned,  its  modes  would  be  sine  and  cosine  waves. 
However,  the  presence  of  mistiming  alters  the  mode  shapes  to 
various  extents.  Some  modes  appear  to  be  nearly  pure 
sinusoids,  while  others  more  heavily  distorted.  Figure  4  shows 
several  representative  sets  of  mode  shape  comparisons  that 
range  from  nearly  tuned-looking  modes  to  modes  that  are  very 
localized.  In  all  cases,  the  modes  from  the  two  methods  agree 
quite  well.  In  addition,  the  natural  frequencies  were  also 
accurately  identified.  Fig.  5.  Thus,  the  traveling  wave  modal 
analysis  method  can  accurately  determine  the  modes  and  natural 
frequencies  of  a  bladed  disk  based  on  its  response  to  a  traveling 
wave  excitation. 

Next,  we  will  demonstrate  that  the  resulting  modes  and 
natural  frequencies  can  be  used  with  FMM  ID  to  identify  the 
mistiming  in  the  bladed  disk.  Since  most  of  the  mistiming  in 
this  fan  was  caused  by  the  attached  masses,  to  a  large  extent 
the  mistuning  is  known.  Therefore,  we  will  use  these  mass 
values  as  a  benchmark  with  which  to  assess  the  accuracy  of  the 
FMM  ID  results. 

Since  we  wish  to  use  the  mass  values  as  a  benchmark,  we 
must  isolate  the  mistuning  caused  by  the  masses  from  the 
inherent  mistuning  in  the  fan.  Therefore,  we  first  performed  a 
standard  SISO  modal  analysis  on  the  rotor  with  the  masses 
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Figure  4:  Comparison  of  the  representative  mode  shape 
extracted  from  the  traveling  wave  response  data  with 
benchmark 
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Figure  5:  Comparison  of  the  natural  frequencies  extracted 
from  the  traveling  wave  response  data  with  benchmark 

removed,  and  used  the  resulting  modal  data  as  input  to  FMM 
ID.  This  resulted  in  an  assessment  of  the  BBR’s  inherent 
mistuning,  expressed  as  a  percent  change  in  each  sector’s 
frequency 

Next,  we  performed  an  FMM  ID  analysis  of  the  modes  and 
frequencies  extracted  from  the  traveling  wave  response  of  the 
rotor  with  mass-mistuning.  The  resulting  mistiming  represents 
the  total  effect  of  the  masses  and  the  IBR’s  inherent  mistuning. 
To  isolate  the  mass  effect,  we  then  subtracted  the  rotor’s 
nominal  mistuning.  Again,  the  resulting  mistuning  was 
expressed  as  a  percent  change  in  each  sector’s  frequency. 

In  order  to  compare  these  mistuning  values  with  the  actual 
masses  placed  on  the  blade  tips,  we  must  first  translate  each 
sector  frequency  change  into  its  corresponding  mass.  A 
calibration  curve  to  relate  these  two  quantities  was  generated 
through  two  independent  methods.  First,  the  calibration  was 
determined  through  a  series  of  finite  element  analyses  in  which 
we  placed  known  mass  elements  on  the  tip  of  a  blade,  and  used 
the  finite  element  model  to  directly  calculate  their  effect  on  the 
corresponding  sector’s  frequency.^  While  this  method  is 
sufficient  in  this  case,  there  are  often  times  when  a  finite 
element  model  is  not  available.  For  such  cases,  a  similar 
calibration  curve  can  be  generated  experimentally  by  varying 
the  mass  on  a  single  blade,  and  repeating  the  FMM  ID  analysis. 
This  experimental  method  was  performed  as  an  independent 
check  of  the  calibration.  Both  approaches  gave  very  similar 
results,  Fig.  6,  For  the  range  of  masses  used  in  this  experiment, 
we  found  that  mass  and  sector  frequency  change  are  linearly 
related.  This  calibration  curve  was  then  used  to  translate  the 


^  A  single  blade  disk  sector  of  the  tuned  bladed  disk  with  cyclic 
symmetric  boundary  conditions  applied  to  the  disk  is  used  in  this  calculation. 
Changing  the  phase  in  the  cyclic  symmetric  boundary  condition  only  had  a 
slight  effect  on  the  results.  The  results  given  in  Figure  5  are  representative  and 
corresponded  to  a  phase  constraint  of  90  degrees. 


6 


Copyright  ©  2003  by  ASME 


identified  sector  frequency  changes  into  their  corresponding 
masses. 

Figure  7  shows  the  comparison  between  the  mass  mistuning 
identified  through  FMM  ID  with  the  values  of  the  actual  masses 
placed  on  each  blade  tip.  The  agreement  is  quite  good.  Thus, 
by  combining  the  traveling  wave  modal  analysis  method  with 
FMM  ID,  we  can  determine  the  mistuning  in  a  bladed  disk  from 
its  traveling  wave  response 


Figure  6:  Calibration  curve  relating  the  effect  of  a  unit 
mass  on  a  sector’s  frequency  deviation 


Blade  Number 

Figure  7:  Comparison  of  the  mistuning  from  the  traveling 
wave  system  identification  method  with  benchmark  value 

3.2  Rotating  Test  Case 

In  the  previous  example,  we  verified  the  traveling  wave 
modal  analysis  method  and  showed  that  the  resulting  modes  and 
natural  frequencies  could  be  used  to  determine  the  bladed  disk’s 
mistuning.  This  was  a  useful  benchmark  experiment  since  the 
rotor  was  excited  while  it  was  stationary,  and  thus  we  could 


make  very  accurate  measurements  using  a  SLDV.  'However,  if 
the  method  is  to  be  applicable  to  conventional  bladed  disks,  we 
must  be  able  to  make  response  measurements  under  rotating 
conditions.  This  second  test  case  is  intended  to  assess  if  the 
measurement  techniques  commonly  used  in  rotating  tests  are 
sufficiently  accurate  to  be  used  with  FMM  ID  to  determine  the 
mistuning  in  a  bladed  disk. 

^  For  this  example,  we  considered  the  fan  shown  in  Fig.  8. 
To  obtain  a  benchmark  measure  of  the  rotor’s  mistuning  in  its 
first  bending  modes,  we  used  an  impact  hammer  and  a  laser 
vibrometer  to  perform  a  SISO  modal  analysis  test.  The 
resulting  modes  and  natural  frequencies  were  then  used  as  input 
to  FMM  ID  to  determine  the  fan’s  mistuning. 


Figure  8:  Fan  studied  in  rotation  tests 


Next,  the  fan  was  tested  in  the  spin  pit  facility  at  NASA 
Glenn  Research  Center.  The  system  uses  an  array  of  permanent 
magnets  to  generate  an  eddy  current  excitation  that  drives  the 
blades.  The  blade  response  is  then  measured  with  an  NSMS 
system.  For  this  test,  ffie  fan  was  driven  with  a  7E  excitation, 
over  a  rotational  speed  range  of  1550  to  1850  RPM.  The  test 
was  performed  twice,  at  two  different  acceleration  rates.  The 
NSMS  signals  were  then  processed  to  obtain  the  amplitude  and 
phase  of  each  blade  as  a  function  of  its  excitation  frequency, 
Fig.  9.  The  NSMS  system  measures  the  amplitude  and  phase  of 
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(a)  Fast  Acceleration 


(b)  Slow  Acceleration 

Figure  9:  Tracking  plot  of  blade  amplitudes  measured  in  spin 

pit 


each  blade  once  per  revolution.  Thus,  the  data  taken  at  the 
slower  acceleration  rate  has  a  higher  frequency  resolution  than 
that  obtained  from  the  faster  acceleration  rate.  In  both  cases, 
the  data  is  significantly  noisier  than  the  measurements 
obtained  in  the  previous  example  using  an  SLDV 

Next,  we  applied  the  traveling  wave  system  identification 
method  to  extract  the  mode  shapes  from  the  response  data. 
First,  we  transformed  the  measurements  to  the  modal  analysis 
domain  by  using  (5),  and  extracted  the  mode  shapes  and 
natural  frequencies  with  MODENT.  The  extracted  modes 
were  then  transformed  back  to  the  physical  domain  through 
equation  (4).  Finally,  the  resulting  modes  and  frequencies 
were  used  as  input  to  FMM  ID  to  identify  the  fan’s  mistuning. 


(a)  Fast  Acceleration 


Figure  10:  Comparison  of  the  mistuning  determined  through 
the  traveling  wave  system  identification  method  with 
benchmark  values 


The  mistuning  identified  from  the  two  spin  pit  tests  was 
then  compared  with  the  benchmark  values,  as  shown  in  Fig. 
10.  In  the  case  of  the  faster  acceleration  rate,  we  were  able  to 
identify  the  trends  of  the  mistuning  pattern,  but  were  not  able 
to  accurately  determine  the  mistuning  values  for  all  blades. 
Fig.  10(a).  The  key  limitation  in  this  analysis  was  our  ability 
to  extract  accurate  mode  shapes  from  data  with  such  coarse 
frequency  resolution.  However,  the  frequency  resolution  of 
the  data  measxired  at  a  slower  acceleration  rate  was  3  times 
higher  than  the  previous  case.  Thus,  when  FMM  ID  was 
applied  to  this  higher  resolution  data  set,  the  agreement 
between  the  traveling  wave  based  ID  and  the  benchmark 
values  was  significantly  improved,  as  shown  in  Fig.  10(b). 
These  results  are  very  encouraging.  They  suggest  that  with 
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adequate  frequency  resolution,  NSMS  measurements  can  be 
used  to  determine  the  mistuning  of  a  bladed  disk  under 
rotating  conditions. 

4.  CONCLUSIONS 

A  method  for  extracting  a  bladed  disk’s  mode  shapes  and 
natural  frequencies  from  its  response  to  a  traveling  wave 
excitation  was  presented.  By  using  the  resulting  information 
as  input  to  the  FMM  ID  system  identification  method,  we  were 
in  turn  able  to  determine  the  mistuning  in  the  bladed  disk. 
Thus,  this  approach  allows  us  to  perform  system  identification 
of  bladed  disks  from  traveling  wave  response  measurements. 
Furthermore,  it  was  demonstrated  that  the  measurement 
technology  that  is  commonly  used  in  rotating  tests,  NSMS, 
may  be  used  for  this  purpose. 

There  are  a  number  of  advantages  to  performing  system 
identification  based  on  a  bladed  disk’s  response  to  a  traveling 
wave  excitation.  First,  it  allows  us  to  use  data  taken  in  a  spin 
pit  or  stage  test  to  determine  a  rotor’s  mistuning.  In  this  way, 
the  mistuning  we  identify  will  include  all  effects  present  during 
the  test  conditions,  i.e.  centrifugal  stiffening,  gas  loading, 
mounting  conditions,  as  well  as  temperature  effects.  In 
particular,  it  will  allow  xis  to  assess  the  effect  that  centrifugal 
loading  of  the  attachments  has  on  the  mistuning  in 
conventional  bladed  disks.  Thus,  this  technique  extends  the 
applicability  of  FMM  ID  from  IBRs  to  conventional  bladed 
disks. 

In  addition,  the  traveling  wave  system  identification 
method  allows  us  to  accurately  assess  the  mistuning  of  some 
bladed  disks  in  a  far  more  efficient  manner.  In  order  to 
determine  a  rotor’s  mistuning,  FMM  ID  uses  its  mode  shapes 
and  natural  frequencies  as  input.  Although  FMM  ID 
theoretically  only  needs  measurements  of  one  or  two  modes, 
the  method’s  robustness  and  accuracy  is  greatly  improved 
when  more  modes  are  included.  For  certain  bladed  disks,  a 
single  traveling  wave  excitation  can  be  used  to  measure  more 
modes  than  would  be  possible  from  a  single  point  excitation 
test.  For  example,  consider  a  highly  mistuned  rotor  that  has  a 
large  number  of  localized  modes.  It  is  often  hard  to  excite  all 
of  these  modes  with  only  one  single  point  excitation  test, 
because  the  excitation  source  will  likely  be  at  a  node  of  many 
of  the  modes.  Therefore,  if  we  wish  to  detect  all  of  the  mode 
shapes,  the  test  must  be  repeated  at  various  excitation  points. 
However,  if  we  were  to  drive  the  system  with  a  traveling  wave 
excitation,  we  can  generally  excite  all  localized  modes  with 
just  a  single  engine  order  excitation.  This  difference  arises 
because  the  more  localized  a  mode  becomes  in  physical 
coordinates,  the  more  extended  it  will  be  in  the  modal  analysis 
coordinate  system.  Thus  in  highly  mistuned  systems,  one 
engine  order  excitation  can  often  provide  more  modal 
information  than  several  single  point  excitations. 


Finally,  the  traveling  wave  system  identification  method 
presented  in  this  paper  could  potentially  form  the  basis  of  an 
engine  health  monitoring  system.  If  a  blade  develops  a  crack, 
its  frequency  will  decrease.  Thus,  by  analyzing  blade 
vibration  in  the  engine,  the  traveling  wave  system 
identification  method  could  detect  a  cracked  blade.  A  health 
monitoring  system  of  this  form  would  use  sensors,  such  as 
NSMS,  to  measure  the  blade  vibration  through  an  accel.  The 
measurements  may  be  filtered  to  isolate  an  engine  order 
response,  and  then  analyzed  using  the  traveling  wave  system 
identification  method  presented  in  this  paper.  This  will  result 
in  a  measure  of  the  rotor’s  mistuning,  which  can  be  compared 
with  previous  measurements  to  identify  if  any  blade’s 
frequency  has  changed  significantly,  thus  identifying  potential 
cracks.  The  key  obstacle  in  automating  this  procedure  is  to 
develop  a  mode  extraction  method  that  does  not  require  user 
interaction.  Although  this  may  not  be  possible  in  general,  it  is 
likely  that  one  could  develop  an  automated  modal  analysis 
method  which  is  tailored  to  a  specific  piece  of  hardware. 
Thus,  the  traveling  wave  system  identification  method  has  the 
potential  to  be  an  extremely  useful  tool  for  engine  health 
monitoring. 

The  focus  of  this  paper  is  on  extracting  the  modes  of  a 
bladed  disk  from  its  response  to  a  traveling  wave  excitation, 
however,  the  general  approach  presented  may  be  extended  to 
any  structure  subjected  to  a  mtilti-point  excitation  in  which  the 
driving  frequencies  are  consistent  from  one  excitation  point  to 
the  next.  A  particular  advantage  to  this  approach  is  that  it 
allows  structures  to  be  tested  in  a  manner  that  more  accurately 
simulates  their  actual  operating  conditions.  Thus,  this 
transformation  technique  for  multi-point  excitation  modal 
analysis  is  not  only  a  useful  tool  for  bladed  disk  system 
identification,  but  potentially  can  be  used  to  address  a  much 
larger  class  of  systems. 
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ABSTRACT 

Various  reduced  order  models  have  been  developed  to  predict  the  forced  response  of  a 
bladed  disk  in  an  engine.  Many  of  these  prediction  codes  rely  on  steady  state 
calculations.  However,  the  high  acceleration  rates  used  in  gas  turbine  engines  can 
produce  transient  behavior  in  the  vibratory  response.  Such  transient  effects  can  be 
significant,  and  may  lead  to  poor  correlation  of  the  predicted  response  with  experimental 
data.  In  addition,  transient  vibratory  response  may  contribute  significantly  to  fatigue 
damage  in  the  engine  and  may  not  be  properly  modeled  using  steady  state  calculations. 

A  new  transient  reduced  order  model  of  bladed  disk  vibration  is  developed.  The  method 
is  an  extension  of  the  Fundamental  Mistuning  Model  (FMM),  and  relies  on  numerical 
evaluation  of  the  differential  equations  of  motion  to  capture  the  transient  effects.  The 
model  is  verified  experimentally. 

The  t  ransient  r  esponse  i  s  f ound  t  o  b  e  f  ar  m  ore  sensitive  t  o  s  mall  e  rrors  i  n  t  he  n  atural 
fi-equencies  of  the  mistuned  system  than  are  the  steady  state  calculations.  As  a  result, 
small  m  odeling  e  rrors  t  hat  p  roduce  n  egligible  e  ffects  i  n  a  s  teady  s  tate  prediction  h  ave 
been  shown  to  produce  20%  amplitude  errors  m  the  transient  solution.  Thus,  accurate 
modeling  is  essential  for  an  effective  transient  forced  response  prediction. 

Since  the  transient  FMM  code  is  far  more  time  consuming  than  the  steady  state  version,  it 
is  important  to  identify  the  acceleration  regimes  in  which  transient  effects  are  significant, 
and  when  they  may  be  n  eglected.  A  series  of  analytical  and  numerical  analyses  were 
performed  to  generate  acceleration  guidelines,  which  may  be  used  to  assess  when  a 
transient  analysis  is  required  and  when  the  simpler  steady  state  analysis  is  sufficient. 


1.  INTRODUCTION 


Bladed  disks  are  generally  designed  to  be  cyclically  symmetric  with  each  blade 
identical.  However,  due  to  imperfections  in  the  manufacturing  process  as  well  as  wear 
during  operation,  there  are  slight  variations  from  blade  to  blade.  These  small  differences 
have  a  large  effect  on  the  vibratory  response  of  the  disk,  causing  some  blades  to  have 
significantly  larger  amplitude  than  would  be  the  case  in  the  nominal  system.  This  is 
referred  to  as  the  mistuning  problem,  and  has  been  studied  extensively,  as  it  contributes 
to  high  cycle  fatigue  and  the  failure  of  bladed  disks  in  service.  Srinivasan  provides  a 
thorough  review  of  this  topic  in  [  1  ] . 

There  have  been  several  reduced  order  models  developed  to  efficiently  predict  the 
behavior  of  mistuned  bladed  disks  [2-5]. 

However,  one  thing  all  these  models  have  in  common  are  that  they  deal  only  with 
the  steady  state  response  of  rotors.  When  hi^  acceleration  rates  are  used  in  engines, 
transient  behavior  is  observed.  An  example  of  this  type  of  behavior  is  shown  in  Figure  1, 
which  is  data  from  a  spin  pit  test  of  an  IBR  measured  at  NASA  Glenn.  In  a  test  such  as 
this,  the  steady  state  response  of  the  IBR  does  not  reflect  its  actual  behavior.  In  order  to 
correlate  data  such  as  this  with  predictions,  a  model  that  accounts  for  transient  behavior  is 
needed.  The  piupose  of  this  research  was  to  develop  and  validate  such  a  transient  model. 

2.  THEORY 

The  transient  simulation  code  is  based  on  the  Fimdamental  Mistuning  Model 
(FMM),  an  existing  steady-state  reduced  order  model  that  is  applicable  to  isolated 
families  of  modes  [5].  An  extension  of  FMM,  called  FMM  ID,  is  a  completely 
experimental  method  of  determining  the  tuned  system  frequencies  and  mistuning  for  a 
rotor  [6,7]. 

When  applying  this  method  to  a  forced  response  calculation,  FMM  ID  is  first  used 
to  determine  the  bladed  disk’s  mistuning  and  tuned  system  frequencies.  These  identified 
parameters  are  used  along  with  FMM  to  calculate  the  disk’s  mistuned  modes.  A  modal 
summation  is  then  used  to  calculate  the  forced  response. 

When  FMM  is  being  applied  to  steady  state  response,  there  is  a  closed  form 
solution  for  the  response  of  each  mode  and  the  modes  can  be  summed  to  determine  the 
response  of  the  system.  In  the  case  of  transient  excitations  the  solution  also  can  be 
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calculated  by  summing  modes,  but  the  modal  response  must  be  calculated  numerically. 
Consider  the  case  when  the  blades  are  excited  harmonically,  but  the  frequency  of  the 
excitation  is  a  function  of  time.  Then,  if  the  modes  are  normalized  so  that  they  have  a 
modal  mass  of  one,  the  modal  equations  of  motion  have  the  form: 

+  (1) 
For  purposes  of  discussion  assxime  that  the  excitation  frequency  is  a  linear  function  of 
time  ,  i.e.  o)(t)  =  cOg  +ft.  When  this  expression  is  substituted  into  (1),  a  second  order 

time  term  appears  in  the  exponential.  As  there  is  no  closed  form  solution  to  this  equation, 
it  must  be  solved  numerically. 


3.  EXPERIMENTAL  TEST  CASES 

This  method  was  validated  by  comparing  it  with  experiments  on  two  separate 
disks,  one  with  a  single  blade,  and  one  with  18  blades.  Photographs  of  each  disk  are 
shown  in  Figures  2  and  3 . 

For  the  single  blade  test  case,  the  vibration  was  excited  acoustically  with  a  two 
second  sine  sweep  through  the  blade’s  first  bending  resonance.  A  laser  vibrometer  was 
used  to  measure  and  record  the  time  history  of  the  blade’s  response.  A  plot  of  this  time 
history  is  shown  in  Figure  4.  Superimposed  on  this  is  a  plot  of  the  response  envelope, 
calculated  by  determining  the  magnitude  of  the  response  at  each  instant.  For  the 
remainder  of  the  plots  in  this  paper,  only  this  response  envelope  will  be  shown. 

It  may  be  noted  from  this  figure  that  there  is  a  beating  phenomenon  present  in  the 
response.  This  occurs  due  to  the  interaction  of  the  mode’s  natural  frequency  and  the 
frequency  of  the  excitation  source.  The  transient  portion  of  the  solution  to  the  equation 
of  motion  always  oscillates  at  the  natural  frequency,  while  the  steady  state  portion  of  the 
solution  oscillates  at  the  excitation  frequency.  Because  these  frequencies  are  close 
together,  a  beating  phenomenon  is  observed.  When  multiple  modes  are  summed  this  type 
of  beating  becomes  very  important  in  establishing  the  overall  system  response.  Since  the 
beating  depends  on  the  spacing  of  the  natural  frequencies,  it  will  be  shown  that  the 
transient  response  is  sensitive  to  variations  in  the  frequencies  of  the  mistuned  modes  in 
quite  a  different  manner  than  is  the  steady  state  response. 
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In  the  first  experiment,  standard  modal  testing  was  used  to  measure  the  natural 
frequency  and  damping  of  the  blade’s  first  bending  mode.  This  information  was  then 
used  as  input  into  the  transient  simulation  code  to  simulate  the  response.  A  comparison 
of  the  m  easured  and  s  imulated  r  esponses  i  s  s  hown  i  n  F  igure  5 .  It  w  as  n  ot  p  ossible  t  o 
directly  measure  the  magnitude  of  the  excitation  force  used  in  the  experiments. 
Consequently,  the  value  used  in  the  simulation  was  chosen  to  minimize  the  error  between 
the  measured  and  simulated  results.  Clearly,  there  is  good  agreement  between  the 
experimental  and  analytical  results  for  the  single  blade  case. 

In  order  to  verify  the  prediction  code  for  a  more  complex  and  realistic  test 
involving  multiple  modes,  we  used  an  eighteen  bladed  disk.  For  this  test,  we  used  a 
magnetic  excitation  source  applied  to  a  single  blade,  and  did  a  two  second  sine  sweep 
through  the  disk’s  first  bending  mode  resonant  frequencies.  Again,  a  laser  vibrometer 
was  used  to  measure  the  response  of  each  blade. 

This  rotor’s  response  was  then  simulated  using  the  transient  code.  FMM  ID  was 
used  to  determine  the  mistuning  and  tuned  frequencies  in  the  IBR  and  then  FMM  was 
used  to  calculate  the  system’s  mistuned  modes  and  frequencies.  Each  mode’s  response 
was  calculated  and  summed  to  get  the  overall  system  response. 

A  comparison  of  the  results  from  a  representative  blade  is  shown  in  Figure  6. 
Again,  the  results  match  well,  showing  that  the  transient  analysis  correctly  predicts  the 
response  of  systems  with  multiple  modes. 

4.  TRANSIENT  SOLUTION  SENSITIVITY 

A  unique  characteristic  of  transient  behavior  is  the  sensitivity  of  the  blade 
response  to  small  variations  in  the  natural  frequencies  of  the  mistuned  system.  To 
demonstrate  this  effect,  we  changed  three  of  the  natural  frequencies  from  the  18  blade 
case  in  our  numerical  solution.  Specifically,  the  frequencies  of  modes  seven  and  nine 
were  increased  by  0.25%,  and  mode  8  was  decreased  by  0.25%.  The  mode  shapes  and 
damping  were  unchanged.  The  steady  state  frequency  response  of  a  representative  blade 
is  shown  in  Figure  7,  with  the  original  and  perturbed  plots  overlaid.  Clearly,  the  change 
in  frequencies  had  a  negligible  effect  on  the  steady  state  response.  In  addition,  the 
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frequency  changes  had  a  relatively  small  effect  on  the  transient  responses  of  individual 
modes,  refer  to  Figure  8.  However,  consider  the  transient  response  of  the  blade  shown  in 
Figure  9.  There  are  three  lines  on  this  plot,  one  from  the  experiment,  one  from  the 
original  transient  calculation,  and  one  calculated  with  the  perturbed  natural  frequencies. 
To  understand  this,  note  from  the  modal  responses  that  each  mode  peaks  at 
approximately  one  second  after  the  exeitation  begins,  i.e.  the  midpoint  on  the  time  scale. 
Until  that  time  the  perturbed  response  follows  the  original  simulation  fairly  well. 
However,  from  one  second  on,  the  perturbed  response  looks  markedly  different.  This  is 
shown  in  Figure  10,  which  expands  the  second  half  of  Figure  9.  The  change  is  due  to  a 
shift  in  the  phase  of  the  beating  of  the  three  perturbed  modes  caused  by  the  small  changes 
in  the  natural  frequencies.  Thus,  when  these  modes  are  summed,  the  phase  changes 
produce  a  large  change  in  the  individual  blade’s  response. 


5.  ACCELERATION  RATES  THAT  CAUSE  LARGE  TRANSIENT  EFFECTS 

Transient  behavior  is  important  when  the  engine’s  acceleration  rate  exceeds  a 
critical  value.  A  dimensionless  acceleration  rate  can  be  defined  in  terms  of  the  natural 
frequency  of  a  mode. 


From  our  simulations  we  have  found  that  transient  effects  become  important  if  a  >a^ 
where 


or. 


(3) 


and  ^  is  the  damping  ratio  for  the  mode.  This  value  is  comparable  to  the  critical 
acceleration  rate  used  for  frequency  scans  in  modal  testing  [8,9].  A  representative  value 
of  a  fast  acceleration  rate  used  in  a  military  engine  is  6,000  rpm  in  4  seconds.  If  you 
assume  a  natural  frequency  of  500  Hz.  and  a  third  engine  order  excitation,  this  works  out 
to  a  dimensionless  acceleration  rate  of  around  5x10"^.  Assuming  a  damping  ratio  of 
0.2%,  the  critical  acceleration  rate  would  be  approximately  1.33  x10”®.  Consequently, 
the  actual  acceleration  rates  that  occur  in  military  engines  are  significantly  larger  than 
that  r  equired  toe  ause  s  ignificant  t  ransient  e  ffects.  T  herefore,  t  ransient  effects  w  ill  b  e 
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important  in  the  response.  To  illustrate  this  point,  a  numerical  simulation  of  the  18 
bladed  disk  was  perfonned  in  which  the  damping  was  changed  to  0.2%  and  the 
dimensionless  acceleration  rate  set  to  6.3x10“^,  in  order  to  match  representative  engine 
conditions.  T  he  r  esponse  o  f  a  r  epresentative  m  ode  i  s  d  epicted  i  n  F  igure  1 1 .  Clearly, 
transient  beating  effects  are  very  prominent  in  the  response. 

6.  CONCLUSIONS 
The  conclusions  are  that: 

•  A  method  of  predicting  the  transient  response  of  bladed  disks  was  developed. 
This  method  uses  a  numerical  approach  to  integrate  the  modal  equation  of 
motions.  As  input  data,  it  uses  the  mistuned  modes  and  natural  frequencies 
determined  by  FMM  ID. 

•  The  approach  was  shown  to  be  accurate  for  both  a  simple  single  mode  system, 
and  integrally  bladed  disk  with  eighteen  blades. 

•  It  was  demonstrated  that  transient  response  is  more  sensitive  to  small  errors  in 
the  system’s  natural  frequencies  than  are  the  standard  steady-state  calculations. 
This  occurs  because  of  interactions  between  multiple  modes  beating  at  slightly 
different  frequencies. 

•  Guidelines  were  developed  to  determine  imder  what  conditions  transient 
calculations  are  required. 
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FIGURES 


Figure  1 :  Transient  results  from  NASA 
spin  pit  test. 


Figure  2:  Single  blade  rotor 


Figure  3:  Eighteen  blade  rotor 


Figure  4;  Time  response  from  single  blade 
test  case,  with  overlaid  response  envelope. 


Figure  5:  Comparison  of  measurement  with 
simulation  for  single  blade 
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(a) 


(b) 


Figure  6:  Comparison  of  measurement  with  simulation  for  two  representative  blades  from 

18  blade  test  case. 


Figure  7:  Comparison  of  actual  and 
perturbed  steady-state  response  for  1 8 
blade  test  case. 


Figure  8:  Comparison  of  measured,  simulated,  and  perturbed  response  for  two  of  the 
perturbed  modes  from  1 8  blade  test  case. 
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response 


Figure  9:  Comparison  of  measured,  Figure  10:  Zoom  in  on  one  second  to  two 

simulated,  and  perturbed  response  for  a  second  region  of  figure  9. 
representative  blade  from  18  blade  test 
case. 


time  -  sec. 


Figure  1 1 :  Transient  behavior  with  realistic 
accel  rate  and  damping  level. 
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EXPERIMENTAL  STUDY  OF  THE  FUNDAMENTAL  MISTUNING  MODEL  (FMM)  FOR 

PROBABILISTIC  ANALYSIS 
M.R.  Rossi,  D.M.  Feiner,  and  LH.  GrifFm 
Department  of  Mechanical  Engineering 
Carnegie  Mellon  University 

ABSTRACT 

FMM  is  a  reduced  order  model  for  efficiently  calculating  the  forced  response  of  a  mistuned  bladed  disk. 
FMM  ID  is  a  companion  program  which  determines  the  mis  tuning  in  a  particular  rotor.  Together,  these 
methods  provide  a  way  to  acquire  data  on  the  mistuning  in  a  population  of  bladed  disks,  and  then  simulate  the 
forced  response  of  the  fleet.  This  process  is  tested  experimentally,  and  the  simulated  results  are  compared  with 
laboratory  measurements  of  a  “fleet”  of  test  rotors.  The  method  is  shown  to  work  quite  well.  It  is  found  that 
accuracy  of  the  results  depends  on  two  factors:  the  quality  of  the  statistical  model  used  to  characterize 
mistuning,  and  how  sensitive  the  system  is  to  errors  in  the  statistical  modeling. 

1.  INTRODUCTION 

Bladed  disks  are  generally  designed  to  be  cyclically  symmetric.  However,  in  practice  manufacturing 
effects,  non-uniform  material  properties,  and  wear  cause  each  blade  to  be  slightly  different  from  the  rest. 
These  blade-to-blade  variations  are  known  as  mistiming.  The  resonant  amplitudes  of  turbine  blades  are  very 
sensitive  to  these  small  variations  in  the  blade  properties.  Therefore,  mistuning  can  significantly  amplify  the 
vibratory  response  of  some  blades,  and  cause  them  to  fail  from  high  cycle  fatigue.  Srinivasan  provides  a 
thorough  review  of  the  topic  [1]. 

One  area  of  mistuning  research  has  focused  on  the  development  reduced  order  models  to  efficiently 
predict  the  forced  response  of  a  mistuned  bladed  disk.  A  variety  of  reduced  order  models  have  been  developed 
by  researchers  at  Carnegie  Mellon  University  [2-3],  the  University  of  Michigan  [4],  and  Imperial  College  [5]. 
Although  these  methods  have  been  shown  to  agree  extremely  well  with  finite  element  simulations  of  a  full 
mistuned  rotor,  some  have  had  difficulty  predicting  the  response  of  actual  hardware  [6].  These  results  suggest 
that  the  source  of  the  error  may  lie  in  our  inability  to  determine  the  correct  input  parameters  to  the  ROMs. 

Therefore,  researchers  subsequently  developed  methods  to  accurately  measure  the  mistuning  in  a  bladed 
disk.  The  first  advanced  mistuning  identification  method  was  created  by  Judge  and  Pierre  [7].  Their  technique 
uses  rneasurements  of  the  bladed  disk  system  as  a  whole  to  infer  the  mistuning  of  individual  blades.  More 
recently,  Kim  and  Griffin  developed  a  similar  technique  which  is  applicable  to  veering  regions  and  high 
frequency  modes  [8],  However,  the  techniques  of  Judge  and  Kim  require  a  finite  element  model  of  the  system, 
and  significant  analysis  to  identify  the  mistuning  in  an  DBR.  A  much  simpler  approach  is  the  FMM  ID  method 
developed  by  Feiner  and  Griffin  [9,10]. 

FMM  ID  is  based  on  the  Fundamental  Mistuning  Model  (FMM),  a  simple  reduced  order  model  for 
mistuned  bladed  disks  [3].  Like  the  methods  of  Judge  of  Kim,  FMM  ID  also  uses  measurements  of  the  whole 
assembly  to  infer  blade  frequencies,  but  it  is  completely  experimental.  FMM  ID  does  not  require  a  finite 
element  model.  This  makes  FMM  ID  easier  to  use  than  the  other  methods.  Furthermore,  we  have  found  that 
in  practice  it  is  extremely  difficult  to  obtain  a  finite  element  model  that  accurately  reflects  a  component’s  true 


geometry  and  boundary  conditions.  Such  inaccuracies  may  lead  to  errors  when  using  the  methods  of  Judge 
and  Kim.  However,  since  FMM  ID  relies  only  on  experimental  data,  it  is  not  hindered  by  the  quality  of  a  finite 
element  model. 

FMM  ED  can  also  experimentally  determine  the  natural  jfrequencies  that  a  bladed  disk  would  have  if  it 
were  tuned.  The  combination  of  tuned  frequencies  and  mistuning  provide  enough  information  for  FMM  to 
predict  the  forced  response  of  a  bladed  disk  based  solely  on  experimentally  measured  data. 

Both  FMM  and  FMM  ID  have  been  verified  experimentally  in  deterministic  calculations  [10].  However, 
it  has  also  been  proposed  that  these  methods  may  be  used  for  probabilistic  analysis.  The  idea  is  to  use  FMM 
ID  t  o  a  cquire  d  ata  o  n  t  he  m  istuning  a  nd  t  uned  s  ystem  f requencies  i  n  a  p  opulation  o  f  b  laded  d  isks.  T  hen, 
perform  Monte  Carlo  simulations  with  FMM  to  assess  the  response  of  the  fleet.  In  this  paper,  we  test  this 
process  experimentally,  and  then  compare  the  simulated  results  with  measurements  of  a  “fleet”  of  test  rotors. 

This  paper  is  organized  as  follows.  In  Section  2,  we  summarize  the  FMM  and  FMM  ID  methods.  Then, 
Section  3  describes  the  benchmark  experiments.  In  Section  4  we  discuss  the  probabilistic  analysis,  and 
compare  our  simulation  results  with  experimental  data.  Finally,  the  key  results  are  summarized  in 
Conclusions. 

2.  FMM  AND  FMM  ID 

2.1  FMM 

FMM  is  a  simple  reduced  order  model  of  mistuned  bladed  disk  vibration.  The  method  is  a  simplification 
of  the  Subset  of  Nominal  Modes  theory  (SNM)  developed  by  Yang  and  Griffin  [2]  and  is  designed  for  use  in 
low  fi'equency  modes  such  as  first  bending  and  first  torsion.  One  of  the  advantages  of  this  approach  is  that  it 
reduces  the  mistuning  problem  to  it  most  basic  elements.  As  a  result,  FMM  requires  a  minimum  number  of 
input  parameters,  and  it’s  extremely  easy  to  use.  This  large  simplification  also  makes  FMM  extremely 
efficient.  When  performing  Monte  Carlo  simulations  of  forced  response,  FMM  can  simulate  the  response  of 
about  200  disks  per  second  on  a  2  GHz  PC. 

The  FMM  method  only  requires  two  sets  of  input  parameters  to  calculate  the  effect  of  mistuning  on  the 
mode  shapes  and  frequencies  of  a  bladed  disk.  Consider  the  eigenvalue  problem  solved  by  FMM  to  calculate 
mistuned  modes  and  frequencies, 

[s:^+a)p,=o,]p^  (1) 

Notice  that  the  equation  has  only  two  input  matrices,  and  A  .  12*"^  is  a  diagonal  matrix  of  the  tuned 

\  ^ 
system  frequencies  squared.  This  term  describes  the  nominal  system.  The  matrix  A  characterizes  the 

mistuning.  A  is  composed  of  the  blade  frequency  deviations,  which  are  defined  as  the  difference  in  each 
blade’s  frequency  from  the  average  value.  Thus,  FMM  shows  that  the  effect  that  mistuning  has  on  a  system  is 
completely  defined  by  only  two  sets  of  parameters:  the  tuned  system  frequencies,  and  the  blade  frequency 
deviations. 
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This  result  has  a  large  implication  to  probabilistic  analysis  because  it  minimizes  the  number  of  random 
variables  which  must  be  accounted  for  when  calculating  the  response  of  mistuned  systems.  Since  misturiing 
may  be  characterized  by  the  frequencies  of  the  blades,  it  is  not  necessary  to  separately  model  the  variations  in 
blade  geometry  and  material  properties.  All  we  need  to  account  for  is  the  net  effect  of  these  variations  on  the 
blade  frequencies. 

Therefore,  it  is  important  to  be  able  to  accurately  measure  the  frequencies  of  individual  blades.  Blade 
frequencies  are  often  difficult  to  measure  directly,  particularly  in  the  case  of  IBRs,  where  the  blades  cannot  be 
removed  for  individual  testing.  But  FMM  provides  a  method  for  blade  frequency  identification:  FMM  ID. 

2.2  FMM  ID 

Recall  that  the  FMM  eigenvalue  problem,  Eq.  (1),  is  used  to  calculate  the  effect  of  mistuning  on  the  mode 
shapes  and  natural  frequencies  of  a  bladed  disk.  The  equation  takes  as  input,  information  on  the  nominal 
system  as  well  as  the  way  it’s  mistuned.  With  this  data,  the  expression  can  be  solved  for  the  mode  shapes  and 
natural  frequencies  of  the  mistuned  bladed  disk.  However,  we  could  alternatively  solve  the  problem  in  reverse. 
Suppose  we  measured  the  modes  and  natural  frequencies  of  a  mistuned  rotor.  We  could  then  formulate  an 
inverse  problem  to  Eq.  (1)  [9].  The  solution  to  this  inverse  problem  is  the  mistuning  of  each  blade,  as  well  as 
the  natural  frequencies  the  system  would  have  if  it  were  tuned.  This  is  the  basis  of  FMM  ID,  and  it  is  shown 
schematically  in  Fig.  1. 

In  practice,  the  mistuned  system  modes  and  frequencies  are  measured  through  standard  modal  testing 
techniques.  This  involves  measuring  a  set  of  transfer  functions,  and  then  extracting  modes  with  modal  curve 
fitting  software.  For  the  purpose  of  FMM  ID,  the  modes  only  need  to  be  measured  at  one  point  per  blade. 

FMM  ID  does  not  require  any  finite  element  data.  Thus,  it  provides  us  with  a  way  to  determine  all  of  the 
key  mistuning  parameters  experimentally. 

2.3  Probabilistic  Application  of  FMM  and  FMM  ID 

By  measuring  multiple  bladed  disks  of  the  same  design,  FMM  ID  can  collect  data  on  mistuning  and  tuned 
system  frequencies,  which  can  later  be  used  for  probabilistic  analysis. 

Once  the  data  is  collected,  we  can  construct  statistical  models  of  the  mistuning  and  tuned  system 
frequencies.  Then,  we  can  use  those  statistical  models  with  FMM  to  perform  Monte  Carlo  simulations  of  the 
fleet. 

The  goals  of  this  research  are  to  apply  this  technique  to  an  academic  rotor,  and  compare  our  simulations 
with  experimental  data.  This  will  allow  us  to  explore  some  of  the  challenges  of  probabilistic  mistuning 
analysis. 

3.  BENCHMARK  EXPERIMENT 

The  experimental  approach  is  to  first  generate  benchmark  data  with  which  to  compare  our  Monte  Carlo 
simulation  results.  This  requires  measuring  the  forced  response  of  multiple  bladed  disks  of  the  same  design. 
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Next,  we  can  apply  the  probabilistic  procedure  discussed  in  the  previous  section  to  simulate  the  response  of  a 
fleet  of  similarly  constructed  disks. 

The  first  step  in  performing  the  benchmark  experiment  was  to  obtain  a  tuned  bladed  disk,  which  we  could 
later  mistune  in  a  controlled  fashion.  Figure  2  shows  the  academic  BBR  used  for  experiment.  To  tune  the 
structure,  we  first  used  FMM  ID  to  measure  the  firequency  of  every  blade.  Since  the  blades  on  this  disk  have  a 
simple  beam-like  geometry,  we  were  then  able  to  use  beam  theory  to  calculate  the  appropriate  length  change 
for  every  blade  to  compensate  for  its  mistiming.  Finally,  we  trimmed  the  blade  lengths  accordingly.  Figure  3 
shows  the  fi-equency  response  function  (FRF)  of  the  rotor  before  and  after  tuning.  Notice  that  prior  to  tuning, 
the  stmcture’s  mistuning  caused  the  repeated  natural  firequencies  to  split,  producing  additional  peaks  of  the 
FRF.  After  timing,  the  splitting  was  eliminated  in  most  of  the  modes.  Thus,  the  disk  was  successfully  tuned. 

Then,  the  disk  was  mounted  in  the  test  fixture  shown  in  Fig.  4.  The  disk  was  mistuned  by  adding  masses 
to  the  blade  tips.  The  masses  were  selected  to  produce  variations  in  the  blade  frequencies  that  were 
approximately  n  ormally  d  istributed  w  ith  a  s  tandard  deviation  e  qual  t  o  2  %  o  f  t  he  nominal  b  lade  frequency. 
Note  that  a  mean  shift  in  the  mistuning  is  mathematically  equivalent  to  a  mean  shift  in  the  tuned  system 
frequencies.  Therefore,  we  defined  the  mean  mistuning  to  be  zero,  and  measured  the  corresponding  tuned 
system  frequencies  through  FMM  ID. 

We  excited  the  disk  with  an  array  of  electro-magnets  positioned  under  the  blade  tips.  Fig.  4.  The  magnets 
produced  an  engine  style  excitation,  while  the  disk  remained  stationary.  The  engine  style  excitation  system  at 
Carnegie  Mellon  is  similar  to  the  one  developed  by  Jones  and  Cross  at  the  Air  Force  Research  Laboratones 
[11].  The  bladed  disk  was  excited  over  an  appropriate  frequency  range  to  simulate  the  effect  of  an  engine 
order  crossing  with  the  first  bending  modes.  The  vibratory  response  of  the  blades  was  measured  at  each  blade 
tip  by  using  a  scanning  laser  vibrometer.  Fig  5a.  Laser  vibrometers  are  ideal  tools  for  mistuning  measurements 
since  they  are  very  accurate,  non-contacting  sensors  which  don’t  alter  the  mistuning  of  the  structure. 

This  measurement  process  was  repeated  with  10  different  mistuning  patterns,  each  drawn  from  the  same 
normal  distribution.  This  effectively  gave  us  measurements  of  10  different  disks  from  the  same  population. 
On  each  test,  we  recorded  the  peak  amplitude  of  every  blade  over  the  frequency  range  of  interest.  Fig  5b. 
Since  each  disk  has  24  blades,  this  produced  a  total  of  240  peak  amplitude  measurements,  which  we  will  later 
use  for  comparison  with  our  simulation  results.  Furthermore,  every  “disk”  was  tested  with  four  different 
engine  order  excitations:  IE,  3E,  6E,  and  9E.  Thus,  we  will  be  able  to  assess  the  accuracy  of  the  probabiUstic 
analysis  method  over  a  wide  variety  of  excitation  conditions. 


4.  PROBABILISTIC  ANALYSIS 

Next,  we  followed  the  probabilistic  analysis  process  outlined  in  Section  2.3.  We  proceeded  as  if  we  knew 
nothing  about  the  way  the  disks  were  mistuned. 

4.1  Single  Disk  Model 

We  performed  a  modal  analysis  on  one  of  the  10  test  “disks.”  Then,  the  measured  modes  were  used  in 
FMM  ID  to  determine  the  stmcture’s  mistuning  and  tuned  system  frequencies.  Fig  6.  This  represents  the 
cmdest  possible  data  for  forming  a  statistical  model  of  the  bladed  disk  parameters.  In  practice,  it  is  advisable 


4 


to  identify  the  parameters  of  multiple  bladed  disks  to  form  a  reliable  model  of  the  random  variables.  However, 
for  the  purpose  of  this  study,  we  would  like  to  assess  the  effect  of  a  crude  statistical  model  on  the  accuracy  of 
the  subsequent  Monte  Carlo  simulations.  Therefore  we  formed  an  approximate  statistical  model  based  on  this 
limited  set  of  data.  With  only  one  measure  of  the  tuned  system  frequencies,  we  have  no  basis  on  which  to 
model  variability.  Therefore,  the  tuned  frequencies  were  treated  as  fixed.  The  blade  frequency  deviations, 
however,  were  modeled  as  a  random  variable.  Figure  7  shows  a  normal  plot  of  the  24  blade  frequency 
deviations  from  this  disk.  Notice  that  the  data  approximately  falls  on  a  straight  line.  This  indicates  that  data  is 
roughly  normal.  Therefore,  the  mistuning  was  modeled  as  being  normally  distributed  with  a  mean  and 
standard  deviation  given  by  the  sample  values  of  0  and  1.52%  respectively.  However,  it  must  me  noted  that 
there  is  substantial  uncertainty  in  these  parameters.  For  instance,  the  95%  confidence  interval  on  the  standard 
deviation  covers  a  range  from  1,18%  to  2.13%.  That’s  nearly  a  factor  of  2  uncertainty  in  the  model  parameter. 

Based  on  this  rough  model,  we  performed  Monte  Carlo  simulations  of  the  bladed  disk  population  using 
FMM.  These  simulations  were  repeated  for  all  four  engine  orders  measured  in  the  benchmark  experiments.  In 
each  case,  we  simulated  the  forced  response  of  1000  bladed  disks.  The  results  are  shown  in  Fig.  8.  Each  plot 
contains  the  CDF  of  all  240  peak  blade  amplitudes  from  the  experiment,  and  a  corresponding  CDF  constructed 
from  the  simulation  results.  For  clarity,  the  plots  are  shown  on  a  normal  probability  scale.  The  agreement  is 
surprisingly  g  ood  c  onsidering  t  hat  t  he  e  xperimental  C  DFs  o  nly  c  ontain  2  40  d  ata  points,  a  nd  a  re  1  ikely  n  ot 
converged  in  the  tails.  Furthermore,  the  statistical  model  used  in  the  simulations  was  inaccurate.  This 
suggests  that  the  response  is  relatively  insensitive  to  errors  in  the  statistical  model.  To  better  understand  this 
behavior,  we  performed  a  sensitivity  analysis. 

4.2  Sensitivity  Analysis 

A  small  change  in  the  standard  deviation  used  in  our  statistical  models  will  produce  a  shift  in  the  simulated 
CDF,  Fig.  9.  In  general,  the  shift  will  not  be  uniform  over  the  full  range  of  the  CDF,  as  shown  in  the  figure. 
Thus,  one  method  for  measuring  sensitivity  is  to  plot  the  change  in  the  CDF  due  to  a  perturbation  in  the 
standard  deviation.  This  analysis  was  performed  about  a  nominal  standard  deviation  of  2%,  and  was  repeated 
for  all  four  engine  orders.  Fig.  10.  Notice  that  the  6  and  9E  cases  are  nearly  zero  across  the  full  range  of 
probability.  Thus,  this  analysis  suggests  that  this  disk’s  response  to  6E  and  9E  excitations  is  very  insensitive  to 
errors  in  the  standard  deviation.  This  is  consistent  with  the  CDFs  of  Fig.  8.  In  particular,  consider  the  9E 
CDF.  Notice  that  despite  a  large  error  in  the  statistical  model  used  to  generate  the  simulated  curve,  it  agrees 
extremely  well  with  the  experimental  data  from  about  10%  to  90%  cumulative  probability.  The  discrepancy 
seen  in  the  tails  is  most  likely  a  result  of  insufficient  experimental  data  to  produce  a  converged  experimental 
CDF  in  those  regions. 

Next,  consider  the  IE  case.  The  sensitivity  plot  indicates  that  the  system’s  response  to  a  IE  excitation  is 
relatively  sensitive  to  errors  in  the  statistical  model  for  most  of  the  probability  range.  However,  in  the  vicinity 
of  85%  cumulative  probability,  the  IE  line  passes  through  zero  on  the  sensitivity  plot,  and  is  therefore  much 
less  sensitive.  Again,  this  is  consistent  with  the  CDF  plot  for  this  case,  which  shows  good  agreement  between 
the  simulation  and  experiment  aroimd  85%,  yet  larger  discrepancies  away  from  that  area.  This  suggests  that 
much  of  the  error  seen  in  the  low  engine  order  simulations  is  due  to  a  poor  statistical  model.  Therefore,  the 
correlation  should  be  improved  if  we  use  a  higher  quality  model  for  the  variation  in  mistuning. 
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4.3  Ten  Disk  Model 


Next  we  formed  a  much  better  statistical  model  of  the  mistuning  by  using  FMM  ID  to  measure  the  blade 
frequency  deviation  in  all  10  test  disks.  Again,  the  data  was  found  to  be  normally  distributed.  In  this  case,  the 
sample  standard  deviation  was  1.92%,  which  is  much  closer  to  the  true  standard  deviation  of  2%.  Since  this 
sample  standard  deviation  is  based  off  of  10  times  as  much  data  as  the  erode  model,  our  uncertainty  in  the 
parameter  has  been  greatly  reduced.  The  95%  confidence  interval  ranges  from  1.77%  to  2.12%. 

The  Monte  Carlo  simulations  were  then  repeated  with  this  improved  statistical  model.  The  resulting  CDFs 
are  shown  in  Fig.  11.  As  expected,  we  see  substantial  improvement  in  the  correlation  of  the  IE  and  3E 
simulations  w  ith  e  xperimental  d  ata.  Furthermore,  s  imulated  C  DFs  f  rom  the6Eand9Ec  ases  a  re  v  irtually 
identical  to  those  generated  from  a  much  cruder  statistical  model.  This  result  confirms  that  the  6  and  9E  cases 
are  insensitive  to  statistical  modeling  errors.  Again,  the  simulations  agree  well  with  the  experimental 
benchmark.  Therefore,  the  FMM  based  probabilistic  analysis  process  may  be  used  to  accurately  determine  the 
statistical  behavior  of  the  fleet. 

4.4  Sensitivity  Dependence  on  Mistuning  Level 

We  found  throughout  Section  4  that  the  accuracy  of  a  probabilistic  mistuning  analysis  depends  on  two 
factors:  the  quality  of  the  statistical  model,  and  sensitivity  ofthe  system’s  response  to  statistical  modeling 
errors.  As  shown  in  Fig.  10,  the  system’s  sensitivity  is  a  function  of  the  engine  order  of  excitation  as  well  as 
the  probability  range  of  interest.  However,  it  should  be  noted  that  the  sensitivity  regime  is  also  governed  by 
the  level  of  mistuning  in  the  system. 

Consider  the  CDFs  shown  in  Fig.  12.  Each  curve  corresponds  to  the  response  of  a  system  with  a  different 
mistuning  standard  deviation,  Oj  and  02  respectively.  In  the  case  of  ctj,  the  99*  percentile  amplitude  is  about 
1.5.  Yet,  the  99*  percentile  amplitude  for  02  is  slightly  higher.  Therefore,  we  can  plot  the  99*  percentile  of 
the  response  as  a  function  of  the  standard  deviation  of  the  mistuning.  Fig  13.  Figure  13  shows  the  99* 
percentile  for  all  four  engine  orders  considered  in  this  study.  Notice  that  in  the  vicinity  of  2%  mistiming,  the 
6E  and  9E  curves  have  a  near-zero  slope.  Thus,  the  99*  percentile  amplitude  in  these  cases  is  insensitive  to 
small  changes  in  the  mistuning  level  (standard  deviation).  This  is  consistent  with  the  sensitivity  plot  shown  in 
Fig.  10.  However,  if  the  mistuning  was  instead  on  the  order  of  0.5%,  then  Fig.  13  indicates  that  the  response 
to  a  6E  or  9E  excitation  would  be  much  more  sensitive  to  changes  in  the  standard  deviation.  Thus,  a  system’s 
sensitivity  to  errors  in  the  statistical  modeling  depends  on  the  level  of  mistuning. 

5.  CONCLUSIONS 

It  was  shown  that  FMM  and  FMM  ID  may  be  used  for  probabilistic  analysis  of  mistimed  bladed  disks.  The 
process  involves  using  FMM  ID  to  collect  data  on  the  mistuning  and  tuned  frequencies  of  a  population  of 
bladed  disks.  This  data  is  then  used  to  construct  statistical  models  of  the  parameters.  Finally,  we  can  use  those 
statistical  models  with  FMM  to  perform  Monte  Carlo  simulations  of  the  fleet  response.  FMM  is  an  ideal 
physical  model  for  Monte  Carlo  simulations  because  it  is  accurate,  simple  to  use,  and  extremely  efficient. 

The  method  was  verified  experimentally  by  comparing  the  results  of  our  Monte  Carlo  simulations  against 
laboratory  measurements  of  mistimed  disks.  The  FMM  approach  worked  very  well.  We  found  that  the 
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accuracy  of  the  method  depends  on  both  the  quality  of  the  statistical  model,  and  the  sensitivity  of  the  system’s 
response  to  errors  in  the  statistical  modeling.  The  sensitivity  regime  may  be  assessed  through  the  sensitivity 
analyses  discussed  in  Sections  4.2  and  4.4.  The  efficiency  of  FMM  makes  these  analyses  fast  and  easy  to 
perform.  If  it  is  found  that  the  system  is  sensitive,  then  the  statistical  models  may  need  to  be  improved  to 
ensure  an  accurate  simulation.  Such  improvements  can  be  made  by  using  FMM  ID  to  measure  the  mistuning 
of  additional  hardware.  Conversely,  additional  testing  may  not  be  necessary  on  systems  that  are  insensitive  to 
modeling  errors. 

FMM  and  FMM  ID  were  experimentally  shown  to  be  effective  tools  for  probabilistic  analysis  of  mistuned 
bladed  disks. 
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Figure  1 :  Schematic  representation  of  the  relation 
between  FMM  and  FMM  ID. 


Figure  2:  TestIBR. 


(b)  After  Tuning 
Figure  3:  FRF’s  of  the  test  IBR  before  and  after  tuning. 
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Figure  5;  Representative  measurements  of  one  disk  configuration,  driven  by  a  IE  excitation. 
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Sector  Frequency  Deviation 
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(b)  Tuned  system  frequencies 


Sector  Frequency  Deviation 

Figure  7:  Normal  plot  of  blade  frequency  deviations  of  one  disk 
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Figure  8:  Comparison  of  the  experimental  and  simulated  CDFs  of  the  peak  blade  amplitudes 


9 


Figure  9:  Change  in  CDF  due  to  perturbation  in 
standard  deviation 
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Figure  1 0:  Sensitivity  of  CDF  to  perturbation  in 
standard  deviation,  centered  about  a=2% 
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Figure  1 1 :  Comparison  of  the  experimental  and  simulated  CDFs  of  the  peak  blade  amplitudes  based  on  both 

the  crude  and  improved  statistical  models. 


Figure  12:  Change  in  99^  percentile  amplitude  due 
to  change  in  mistuning  standard  deviation 
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Figure  13:  99*  percentile  amplitude  as  a  function  of 
mistuning  standard  deviation 
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